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ABSTRACT 


The  factors  which  govern  the  distribution 
of  electrons  and  ions  in  the  earth's  exosphere 
are  discussed.  The  theory  takes  Into  account 
the  effect  of  the  electric  field  which  arises 
from  charge  separation,  the  centrifugal  force 
arising  from  the  rotation  of  the  earth  and  the 
effect  of  the  earth's  gravitational  field.  It 
is  assumed  that  the  charged  particles  are  con¬ 
strained  to  move  only  along  the  direction  of 
the  earth's  magnetic  lines  of  force.  The 
modifications  that  result  in  the  electron  and 
ion  distributions  when  a  temperature  variation 
is  assumed  along  a  line  of  force  are  also 
considered. 

The  results  predicted  by  the  theory  are 
compared  with  actual  experimental  observations 
of  the  exospheric  plasma  which  have  been  obtained 
in  recent  years  using  whistlers  and  using  top¬ 
side  lonograms  made  by  the  Alouette  satellite. 


ill 


SEL-63-110 


nmirntriiimc 


\ji'*  w 


Page 


Abstract  iii 

Contents  Iv 

Tables  v 

Figures  v 

N  omenc la  ture  xi 

I,  Introduction  1 

II.  The  Problem  and  the  Asstmiptlons  2 

III.  Theoretical  Formulation  of  the  Solution  4 

A,  Force  on  Unit  Mass  4 

B,  Derivation  of  Equilibrium  Distributions  7 

C,  Form  of  Solution  for  Equilibrium  Ion  and 

Electron  Densities  17 

D,  Calculation  of  Temperature-Modified 

Geopotential  Height  z  23 

E,  Calculation  of  Ratio  of  Electron  Densities 

at  Conjugate  Points  26 

1.  Temperature  Constant  in  Each  Hemisphere  26 

2.  Temperature  Changing  Along  a  Line  of 

Force  29 

IV.  Results  of  the  Theory  32 

A,  Introduction  32 

B.  The  Electron  and  Ion  Distributions  33 

V.  Comparison  with  Other  Work  and  Summary  of 

Results  42 

A.  Discussion  42 

B.  Conclusions  49 

Acknowledgements  51 

Appendix  A.  Geometry  of  the  Dipole  Field  52 

Appendix  B.  Note  on  the  Ion  Distribution  Curves  34 

References  58 

Table  1  6l 

Table  2  62 

Figures  8-22  (others  in  text)  63 


SEL-63-110 


Iv 


rp  A  m  pc; 


Table 

Table 

NOTE: 
Fig.  1 

Fig.  2 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 


Page 

1.  Temperature  and  Composition  Assumptions  6l 

2.  Summary  of  Presentation  of  Results  62 


FIGURES 

Figures  1-7  sre  Included  in  the  text. 

Figures  8-22  are  presented  at  the  back  of  the  report. 

.  Gravitational  and  centrifugal  forces  and  fc)  4 
acting  on  unit  mass  located  at  the  point  A 
defined  by  one  of  the  coordinate  pairs  (r,e), 

(6^,0)  or  (9q,s). 

.  The  relative  m.agnitudes  of  the  gravitational  and  6 
centrifugal  forces  at  different  points  along 
a  line  of  force. 

.  To  Illustrate  an  element  of  volume  taken  along  8 
a  line  of  force.  The  partial  pressures  of  the  ions 
and  electrons,  p^  and  Pg,and  the  electric 
field  E, (Mange  I960)  are  also  shown. 

Distribution  of  electron  and  ion  densities  with  21 
the  temperature -modified  geopotential  height,  z, 
for  the  case  C  =  1  and  constant  temperature 
T  =  1000° K. 

.  Variation  of  temperature -modi fled  geopotential  25 
height,  z,  with  altitude,  h,  for  an  isothermal 
exosphere. 

i.  Parameters  used  to  calculate  Reo  st  the  27 

reference  level  at  two  magnetically  conjugate 
points  when  the  temperature  in  the  northern 
and  southern  hemispheres  are  different  but 
constant  in  each  hemisphere. 

’.  Parameters  used  to  calculate  Reo  at:  the  29 

reference  level  at  two  magnetically  conjugate 
points  when  the  temperature  is  taken  as  a 
continuous  function  Tjj(s)  in  the  northern 
hemisphere,  and  T3(s)  in  the  southern 
hemisphere.  At  points  1  and  2,  the  temperatures 
are  the  same. 
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NOTE 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 


Figs.  8-22  are  presented  at  the  end  of  the  report. 

8.  Variation  of  the  ion  composition  at  500  km  with 
the  temperature  at  that  level  (Bauer  1963). 

9a.  Relative  electron  and  ion  densities  as  a  function 
of  the  geopotential  height,  z,  for  an  isothermal 
exosphere.  The  form  of  the  distributions  depends 
strongly  on  the  assumed  base  level  composition  and 
the  assumed  temperature  as  may  be  seen  by  comparing 
the  above  with  the  subsequent  figures. 

9b.  Relative  electron  and  ion  densities  as  a  function 
of  the  geopotential  height,  z,  for  an  Isothermal 
exosphere.  Comparison  with  Fig.  9c  indicates 
the  effect  on  the  distributions  of  changing  the 
temperature  only. 

9c.  Relative  electron  and  ion  densities  as  a  function 
of  the  geopotential  height,  z,  for  an  Isothermal 
exosphere.  Comparison  with  Fig.  9b  indicates 
the  effect  on  the  distributions  of  changing  the 
temperature  only.  Comparison  with  Fig.  9d  reveals 
the  effect  of  changing  the  base  composition  only. 

9d.  Relative  electron  and  ion  densities  as  a  function 
of  the  geopotential  height,  z,  for  an  isothermal 
exosphere.  Comparison  with  Fig.  9c  reveals  the 
effect  of  changing  the  base  composition  only. 

9e.  Relative  electron  and  ion  densities  as  a  function 
of  the  geopotential  height,  z,  for  an  isothermal 
exosphere  with  a  relatively  high  temperature. 

10.  Summary  of  the  results  of  Figs.  9a  -  e  for  the 
relative  electron  density  distributions  in  the 
exosphere  over  a  wide  range  of  conditions. 

11.  The  variation  with  temperature  of  the  levels 

at  which  the  ions  indicated  have  equal  abundances. 
The  broken  lines  correspond  to  the  cases  when 
the  composition  at  the  base  level  depends  on 
the  temperature  there  (through  equations  (65)). 

The  continuous  lines  correspond  to  the  cases  when 
the  composition  at  the  base  level  is  independent 
of  temperature  (Composition  1).  These  results 
are  in  agreement  with  those  pres«.  nted  earlier  by 
Bauer  (I962). 
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FicCi  12b  i  *3^0  illu5trst6i  for  3n  isot-h6rni8l  6Xoaph®r®; 

the  relative  elec cron  density  distributions 
along  a  field  line  with  feet  (at  500  to  above 
the  earth)  at  geomagnetic  latitude  45°  ( see 
Tables  1  and  2). 

Fig.  12b.  To  illustrate,  for  an  isothermal  exosphere, 
the  relative  electron  density  distributions 
along  a  field  line  with  feet  (at  500  km  above 
the  earth)  at  geomagnetic  latitude  55°  (  see 
Tables  1  and  2).  The  broken  lines  indicate 
the  distributions  which  would  have  been  obtained 
if  the  centrifugal  force  had  been  neglected. 

Fig.  12c.  To  Illustrate,  for  an  isothermal  exosphere, 
the  relative  electron  density  distributions 
along  a  field  line  w5th  feet  (at  500  km  above 
the  earth)  at  geomagnetic  latitude  65°  (  see 
Tables  1  and  2).  The  broken  lines  Indicate 
the  distributions  which  would  have  been  obtained 
if  the  centrifugal  force  had  been  neglected. 

Fig.  13.  Schematic  diagram  to  Illustrate  the  variation 

of  temperature  along  a  line  of  force.  Ti(s)  is 
used  for  the  temoerature  distributions  in  the 
northern  ^ winter^  hemisphere  and  T2(s)  for  the 
southern  (summer)  hemisphere. 

Fig.  I4a.  To  illustrate,  for  a  non -Isothermal  exosphere 
(see  Fig.  13)»  the  relative  electron  density 
distribution  along  a  field  line  with  feet  (at 
500  km  above  the  ea-'th)  at  geomagnetic  latitude 
35°  (  see  Tables  '  .nd  2). 

Fig.  l4b.  To  illustrate,  for  a  non-isotherraal  exosphere 
(  see  Fig.  13),  the  relative  electron  density 
distribution  along  a  field  line  with  feet  (at 
500  to  above  the  earth)  at  geomagnetic  latitude 
45°  (  see  Tables  1  and  2). 

Fig.  I4c.  To  Illustrate,  for  a  non -isothermal  exosphere 
(  see  Fig.  I3) ,  the  relative  electron  density 
distribution  along  a  field  line  with  feet  (at 
500  to  above  the  earth)  at  geomagnetic  latitude 
55°  (  see  Tables  1  and  2) .  The  broken  lines 
indicate  the  distributions  which  would  have  been 
obtained  if  the  centrifugal  force  had  been 
neglected. 
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Fig.  I4d.  To  illustrate,  tor  a  non-iso^hermal  exosphere 
(see  Fig.  13) j  the  relative  electron  density 
distribution  along  a  field  line  with  feet  (at 
500  km  above  the  earth)  at  geomagnetic  latitude 
65°  ( see  Tables  1  and  2) .  The  broken  lines 
indicate  the  distributions  which  would  have  been 
obtained  if  the  centrifugal  force  had  been 
neglected. 


Fig.  15a.  The  predicted  ratio  of  the  electron  density  in 

winter  to  that  at  the  same  location  and  altitude 
in  summer  for  a  number  of  different  altitudes. 
Composition  1  (see  Table  l)  is  used  at  ^00  km, 
as  was  done  previously.  The  assumed  temperature 
distribution  is  illustrated  in  Fig.  13* 


Fig.  15b.  The  predicted  ratio  of  the  electron  density  in 

winter  to  that  at  the  same  location  and  altitude 
in  summer  for  a  number  of  different  altitudes. 
Compositions  1  and  2f  (see  Table  l)  are  used  at 
500  km,  as  was  done  previously.  The  assumed 
temperature  distribution  is  illustrated  in  Fig.  13. 

Fig.  16a.  Equatorial  N(h)  profile.  Thp  electron  density  at 

1000  km  was  assumed  to  be  lO^^/cc  (see  Tables  1  and  2). 
The  temperature  distributions  along  a  field  line  are 
given  by  equation  (66)  and  are  illustrated  schemati¬ 
cally  in  Fig.  13. 

Fig.  16b.  N(h)  profile  at  geomagnetic  latitude  30°,  The 

electron  density  at  1000  km  was  assumed  to  be  10^/cc 
(see  Tables  1  and  2).  The  temperature  distributions 
along  a  field  line  are  given  by  equation  (66)  and 
are  illustrated  schematically  in  Fig.  I3. 

Fig.  16c.  N(h)  profile  at  geomagnetic  latitude  60° .  The  electron 
density  at  1000  km  was  assumed  to  be  10^/cc  (see 
Tables  1  and  2).  The  temperature  distributions  along 
a  field  line  are  given  by  equation  (66)  and  are 
illustrated  schematically  in  Fig.  I3. 


Fig.  17.  The  average  quiet  day  electron  density  at  1000  km 
(based  on  the  results  of  Thomas  and  Sader  I963) . 
The  continuous  lines  are  observed  values  and  the 
broken  lines  are  extrapolations  which  are  roughly 
consistent  with  equatorial  Alouette  data. 
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Fig.  l8a .  Theoretically  predicted  equatorial  N(h)  profiles 
based  on  observed  Alouette  data  for  the  electron 
density  at  1000  km  for  summer  days .  Curves  for  a 
wide  variety  of  isothermal  exospheric  temperatures 
and  compositions  are  illustrated  (see  Tables  1  and  2). 

Fig.  l8b.  Theoretically  predicted  equatorial  N(h)  profiles 
based  on  observed  Alouette  data  for  the  electron 
density  at  1000  km  for  summer  nights .  Curves  for 
a  wide  variety  of  Isothermal  e.xcspheric  tempera¬ 
tures  and  compositions  are  illustrated  (see 
Tables  1  and  2) .  Mote  the  marked  increase  in  the 
rate  of  ^all  off  of  electron  density  with  height 
at  great  distances  from  the  earth. 

Fig.  l8c.  Theoretically  predicted  equatorial  N(h)  profiles 
based  on  observed  Alouette  data  for  the  electron 
density  at  lOCO  km  for  winter  days .  Curves  for  a 
wide  variety  of  isothermal  exospheric  temperatures 
and  compositions  are  illustrated  (see  Tables  1  and 
2). 

Fig.  l8d.  Theoretically  predicted  equatorial  N(h)  profiles 
based  on  observed  Alouette  data  for  the  electron 
density  at  1000  km  for  winter  n ight s .  Curves  for  a 
wide  variety  of  isothermal  exospheric  temperatures 
and  compositions  are  illustrated  (see  Tables  1 
and  2) . 

Fig.  l3e.  Summary  of  the  results  of  Figs.  I8a  -  d.  The 

curves  shown  are  the  predicted  theoretical  N(h) 
distribut ions . 

Fig.  193.  Theoretical  N(h)  profiles  predicted  by  a  number  of 
workers.  The  temperatures  foj*  curves  a  -  e  are 
2000,  1000,  750,  1500.  1500  K  respectively. 

Composition  1  (Table  1)  has  been  used  for  curves 
c,  d  and  e.  Curves  c  and  d  are  based  on  Alouette 
data  for  summer  nights  and  curve  e  for  Alouette 
data  for  summer  days.  Note  that  the  slopes  of 
curves  c,  d  and  e  are  substantially  different  from 
the  others  near  1000  km.  Also,  the  slopes  of  curves 
c  and  d  at  great  distances  from  the  earth  are 
different  from  the  others. 

Fig.  19b.  Examples  of  experimental  obseiTyatlons  of  exospheric 
electron  density  (equatorial  profiles). 
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Pig.  20.  Comparison  of  experimental  and  theoretical 

exospheric  N(h)  profiles  in  the  equatorial  plane. 

The  curves  a  -  g  are  those  described  in  Fig.  19a. 

The  diagonally  shaded  area  defines  approximately 
the  region  within  which  the  experimental  profiles 
of  Fig.  19b  are  observed  to  lie.  The  vertical 
shading  near  1000  km  Indicates  the  approximate 
region  in  which  the  Alouette  observations  for 
1000  km  lie. 

Fig.  21.  Comparison  of  theory  and  experiment.  The  theoi’et ical 
curve  was  computed  assuming  a  constant  exospheric 
temperature  of  1000°K  and  ionic  Composition  i  (Table 
1)  at  the  base  level.  The  electron  density  at 
1000  km  was  taken  to  be  that  given  by  Alouette  for 
summer  night  conditions  (Fig.  17) . 

Fig.  22.  Observed  and  predicted  ratios  of  the  electron 
density  at  1000  km  at  conjugate  points  in  the 
winter  and  summer  hemispheres.  The  middle  line 
giving  the  observed  average  quiet  day  ratios  is 
based  on  summer  and  winter  Alouette  observations 
at  Stanford.  The  ratio  is  approximately  0.5  over 
the  latitude  range  considered.  The  theoretically 
predicted  ratios  are  in  reasonably  good  agreement. 
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NOMENCLATURE 


NOTE;  The  suffixes  e  and  i  are  used  to  denote  electrons 
and  Ions  respectively,  and  the  suffix  o  is  used  to 
denote  values  measured  at  the  base  level  usually 
taken  to  be  500  kilometers  above  the  earth, 

■n  electron  density  at  500  km  relative  to  O”^ 

density  at  that  level 

ion  density  at  500  km  relative  to  O'*"  density 
at  that  level 


h 

m 

e 

m. 

i 

m . 


ri  or  N 
e 


n 


eo 


n : 


altitude 

scale  height  of  non-ionized  ion,  at  500  km 
electron  mas=; 
ion  mass 

temperature-weighted  ion  mass  average 
electron  density 

electron  density  at  the  base  level 

ion  density;  1  =  1,2,3,  for  O"*",  He"^  and 
respectively 


Ion  density  at  the  base  level 

P _  ratio  of  electron  densities  at  500  km  above  the 

earth  at  conjugate  points  in  the  two  hemispheres 

,9)  the  distance  measured  along  a  line  of  force  from 
its  foot  (at  500  km)  to  the  point  (0  ,9)  on  the 
field  line  ° 


5’(9  ,e')  the  distance  measured  along  a  line  of  force  from 
its  foot  (at  500  km)  to  the  point  (9  ,9’).  It  is 
used  as  an  upper  limit  of  integration 

the  distance  measured  along  a  field  line  from  the 
‘  geomagnetic  equatorial  plane  to  a  point  on  the  base 

level  at  latitude  9_ 
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T 

^eo 

Tlo 

To 

T(s) 


electron  or  ion  temperatures  when  they  are  assumed 
the  same  in  an  Isothermal  exosphere 

electron  temperature 

electron  temperature  at  500  km 

ion  temperature 

ion  temperature  at  500  km 

electron  or  ion  temperature  at  500  km,  when  they 
are  a8s\imed  equal 

electron  or  ion  temperatures,  when  they  are 
assumed  equal 

geomagnetic  latitude  of  a  point  on  a  field  line 


0  geomagnetic  latitude  of  a  field  line  at  500  km 

°  above  the  earth’s  surface 


0’  geomagnetic  latitude  of  a  point  on  a  field  line, 

when  this  point  is  the  upper  limit  of  an  integration 

coordinates  of  a  point  on  a  magnetic  field  line  in 
terms  of  the  latitude,  0,  of  the  point  and  the 
latitude  0  ,  of  the  foot  (at  500  km)  of  the  particu¬ 
lar  field  °llne  concerned 

z(0  ,0*)  temperature  modified  geopotential  height,  corres¬ 
ponding  to  the  point  [B  ,Q*),  For  an  isothermal 
exosphere,  z  is  sometimes  referred  to  as  the 
geopotential  height 
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I.  INTRODUCTION 


A  great  deal  of  experimental  data  on  the  distribution  of 
ions  and  electrons  in  the  upper  parts  of  the  earth's  ionosphere 
are  becoming  increasingly  available  as  the  new  techniques  re¬ 
quired  to  investigate  the  physical  properties  of  the  exosphere 
at  tnese  heights  are  developed.  An  excellent  brief  report  on 
activities  in  this  field  has  been  given  by  Bordeau  (1963)  and 
the  main  references  to  experimental  work  are  quoted  in  that 
paper. 

In  this  report,  an  account  is  given  of  theoretical  con¬ 
siderations  which  govern  the  distribution  of  electrons  and 
ions  above  500  km  in  this  region  of  space.  These  distributions 
are  derived  on  the  assumptions  that  the  charged  particles  are 
constrained  to  move  only  along  the  direction  of  the  earth's 
magnetic  lines  of  force  under  the  action  of  the  earth's  grav¬ 
itational  field  and  of  the  centrifugal  force  arising  from  the 
rotation  of  the  earth.  The  effect  of  the  electric  field,  E, 
which  arises  from  the  charge  separation  that  occurs  [Mange 
(i960)]  resulting  from  the  tendency  of  electrons  to  rise  with 
respect  to  the  heavier  positive  ions  is  also  taken  into  account. 
The  modifications  that  result  in  the  electron  and  Ion  distri¬ 
butions  when  a  temperature  variation  along  a  line  of  force  is 
assumed, such  that  there  are  relatively  large  temperature  dif¬ 
ferences  between  the  northern  and  southern  hemispheres  in  the 
vicinity  of  the  peak  of  the  F  region, are  also  considered. 

Special  attention  is  given  to  the  distributions  that  result 
when  computed  along  a  line  of  force,  and  the  resulting  differ¬ 
ences  in  vertical  profiles  observed  at  different  latitudes  are 
discussed.  It  is  found  useful  and  convenient  to  present  the 
predicted  ratios  of  electron  density  at  a  given  level  in  the 
northern  hemisphere  to  that  obtained  at  the  corresponding  con¬ 
jugate  point  in  the  southern  hemisphere.  The  theory  predicts 
that,  under  certain  circumstances,  this  ratio,  when  plotted 
against  geomagnetic  latitude,  exhibits  two  maxima,  one  in  very 
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low  latitude?  and  the  other  at  high  latitudes.  The  relative 
magnitudes  and  positions  of  these  maxima  depend  to  a  certain 
extent  on  the  precise  assumptions  that  are  made  about  the 
relative  densities  of  the  constituents  at  the  base  of  the 
exosphere  and  on  the  assumed  temperature  distributions. 

In  Section  III,  the  solution  of  the  problem  outlined  in 
Section  II  is  presented.  In  Section  IV,  the  results  of  the 
theory  to  be  expected  on  certain  reasonable  assumptions  about 
the  physical  properties  of  the  exosphere  are  given.  Also  in 
this  section,  the  predictions  of  the  theory  are  compared  with 
the  results  of  a  number  of  other  workers  who  have  considered 
this  problem  from  a  theoretical  point  of  view  and  the  main 
differences  are  outlined.  In  particular,  reference  is  made 
to  the  calculations  of  Johnson  (i960).  Bates  and  Patterson 
(1961),  Hanson  (1962),  Bauer  (196i,  1963),  Rothwell  (1962), 
Gliddon  (1963)^  and  Hanson  and  Patterson  (1963).  Also  in 
Section  rv,  the  results  predicted  by  the  theory  presented  in 
the  report  are  compared  with  actual  experimental  observations 
which  have  been  maae  in  recent  years.  In  particular,  a 
comparison  is  made  with  the  results  presented  by  workers  who 
have  considered  data  obtained  from  the  Alouette  satellite 
(see,  for  example.  King  (1963)>  Thomas  aid  Sader  (1963)  and 
references  quoted  therein).  The  main  conclusions  are  described 
and  summarized  in  the  last  section. 

II.  THE  PROBLEM  AND  THE  ASSUMPTIONS 

The  electrons  and  ions  in  the  earth’s  exosphere  are  acted 
upon  by  a  number  of  important  forces.  These  include  the  gravi¬ 
tational  attraction  toward  the  earth,  the  centrifugal  force 
due  to  the  rotation  of  the  earth  about  its  geographic  axis, 
and  the  force  due  to  the  electric  field,  E,  resi  Iting  from 
the  tendency  of  electrons  to  rise  with  respect  to  the 
heavier  positive  ions  (Mange  i960).  In  this  region  of 
space,  the  particles  are  in  diffusive  equilibrium  and  obey 
Dalton’s  law  of  partial  pressures  so  that  each  constitjeiit  has 
a  partial  pressure  corresponding  to  that  which  it  would  have  if 
it  were  the  only  constituent  present.  The  distribution  of 
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electrons  and  ions  which  would  arise  under  the  circumstances 
outlined  immediately  above  will  be  modified  if  there  is  a  gra¬ 
dient  of  temperature  along  a  line  of  force.  The  resulting 
effects  are  discussed  in  detail  below.  The  main  assumptions 
are  as  follows: 

1.  The  upper  atmosphere,  above  500  km,  consists  of  neutral 
particles  together  with  a  neutral  mixture  of  singly 
charged  positive  tons  and  electrons  only,  and  these  are^ 
in  diffusive  equilibrium.  The  positive  ions  are  0"^,  He 
and  K"^. 

2.  The  partial  pressure  for  each  species  is  balanced  by  the 
earth’s  gravitational  and  centrifugal  forces  and  the 
electric  field  arising  from  charge  separation. 

3.  The  charged  particles  are  constrained  to  move  only  along 
the  lines  of  force  so  that  the  distributions  along  dif¬ 
ferent  lines  of  force  are  quite  independent  of  each  oth¬ 
er. 

4.  No  electrons  are  produced  by  the  action  of  the  sun’s 
ionizing  radiations  above  500  tai. 

5.  The  rate  at  which  electrons  recombine  Is  so  small  in 
comparison  with  other  effects  that  the  loss  of  electrons 
from  this  cause  can  be  neglected  in  the  calculations. 

6.  The  axis  of  rotation  of  the  earth  coincides  with  the 
magnetic  dipole  axis. 

7.  A  difference  of  temperature  can  exist  between  the  north¬ 
ern  and  southern  hemispheres.  The  ratio  Iq/Hjh  of  the 
temperatures  in  the  southern  and  northern  hemispheres 
respectively  is  not  greater  than  1.5>  consistent  with 
satellite  drag  observations  [King-Hele  and  Walker  (1960)3 

8.  The  re native  densities  of  0^,  He^  and  are  known  at 
the  reference  level  at  5OO  km. 

The  calculations  described  in  the  subsequent  paragraphs 
are  aimed  at  deriving  the  resulting  electron  distributions  a- 
long  a  line  of  force. 
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III.  THEORETICAL  FORMULATION  OF  THE  SOLUTION 
A.  FORCE  ON  UNIT  MASS 

The  total  force  g  acting  on  a  unit  mass  at  a  point  A, 
represented  by  the  coordinates  (r,  6),  (e^,  9)  or  (0^,  s) 
(Fig.  l).  Is  given  by 


f 

G.  1.  GRAVITATIONAL  /LVD  CENTRIFUGAL  FORCES  (7  and  f  )  ACTING  ON  UNIT  MASS 
LOCATED  AT  THE  POINT  A  DEFINED  BY  ONE  OF  THE*C00RDIRATE  PAIRS  (r.9).  (0^,0) 
OR  (e^.s) . 
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where  is  the  total  force  on  a  unit  mass  at  the  reference 

level  at  the  poles,  fi  is  the  angular  velocity  of  rota¬ 

tion  of  the  earth  about  the  geographic  axis,  f  and  f  are 
the  gravitational  and  centrifugal  forces  respectively,  u^ 
and  u  are  unit  vectors  In  the  directions  of  r  and  x,  re- 
;-pectlvely,  and  the  distance  s  Is  measured  along  a  line  of 
force  (Fig.  l).  It  is  clear  that  the  force  on  the  unit  mass 
when  It  Is  over  the  poles  acts  toward  the  center  of  the  earth, 
and  that  at  some  particular  points  which  may  be  determined  on 
each_^lde  of  the  earth  In  the  equatorial  plane,  the  forces  fg 
and  f(,  become  equal  and  opposite.  In  fact.  It  Is  possible 
to  define  a  curve  where  the  magnitudes  of  the  two  forces  (but 
not  their  directions,  except  on  the  equatorial  plane)  are  e- 
qual. 

Figure  2  shows  the  locus  of  points  for  which  fg  and  f^ 
are  equal  In  magnitude  (curves  AB).  The  magnitudes  of  fg 
and  f  and  the  resultant  force  g(r,0)  on  a  unit  mass  are 
drawn  to  ::cale  In  Fig.  2  for  different  positions  along  a  line 
of  force  through  70  degrees  geomagnetic  latitude.  If  it  Is 
assumed  that  charged  particles  can  move  only  along  the  lines 
of  force  of  the  earth's  magnetic  field,  then  it  becomes  impor¬ 
tant  to  consider  the  way  In  which  the  resultant  force  f(r,0) 
obtained  by  considering  only  the  resolved  component  of  g(r,e) 
along  the  direction  of  the  earth's  lines  of  force  varies  in 
space  (Fig.  2).  The  loci  of  points  at  which  f(r, 0)  is  zero 
are  shown  as  the  broken  curves  (CD)  in  Pig.  2.  Referring  to 
Fig.  1  and  Eq.  (i).  It  Is  clear  that  the  algebraic  value  of 
this  force  f  (written  without  the  arrow)  Is  given  by 

r  2 

f(”>®)  =  So — 2  ®  ^  ^  ®  ^  (2) 


From. Fig.  1,  it  Is  seen  that 


cos 


^  sin  Of  = 


tan  a 


2tan  9 


>/  1+tan^a  A+^itan^e 


2  sin  0 

COS0  - - - 

vl+4tan‘^0 


(3s) 
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FIG.  2.  THE  RELATIVE  MAGNITUDES  OF  THE  GRAVITATIONAL  AND  CENTRIFUGAL  FORCES  AT 
DIFFERENT  POINTS  ALONG  A  LINE  OF  FORCE. 


cos  Y  =  sln(a+e)  =  sin6  cosa  +  sinv  cos  S 


=  cosa(sin6  +  tana  cos  6) 

sin9-ftanacos9  ^  sin9+2tan9cos9  _  3  ^in  9 

+  tan^a  ^  tan^9  yi+^tan^"  9 


(3b) 


where  a  is  the  magnetic  dip  at  (r,9),  given  by 

tan  a  =  2  tan  9 


Substituting  the  above  values 
(r,9)  to  (©o^e)  where  9^  is 
reference  level  r^,  and  using 

(£SoCOs'*e 

U  i  p  ^ 

Vl+4tan  9 


in  (2),  changing  the  variables 
the  geomagnetic  latitude  at  the 
(Al)  of  Appendix  A: 

2 

30  r^  , 

- i- - ^ -  COS^9)  (^) 

cos^9  cos  9^ 


Equation  (4)  gives  the  resultant  force  resolved  along  the  di¬ 
rection  of  the  earth’s  magnetic  field  line,  arising  from  the 
combined  action  of  the  earth’s  centrifugal  and  gravitational 
forces.  It  is  clearly  a  function  of  the  geomagnetic  latitude 

9  at  vfnich  the  line  of  force  under  consideration  crosses  the 
o 

reference  level,  and  of  the  distance  from  the  earth’s  surface 
to  the  point  under  consideration  as  represented  by  the  coordi¬ 
nate  6,  Fig.  1.  As  will  be  seen  later,  the  precise  way  in 
which  this  force  varies  from  point  to  point  in  space  is  ex¬ 
tremely  important  for  the  physics  of  the  exosphere. 


B.  DERIVATION  OF  EQUILIBRIUM  DISTRIBUTIONS 

There  is  strong  evidence  both  on  theoretical  grounds 
[Nicolet  (1961)]  and  on  experimental  grounds  [e.g.  Hanson 
(1962)]  for  supposing  that  the  main  ionized  constituents  of 
the  upper  atmosphere  above  the  peak  of  the  P2  layer  are  elec 
trons,  together  with  oxygen,  helium  and  hydrogen  ions  in  the 
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atomic  stato.  In  the  analysis  described  belcw,  it  is  assumed 

that  these  are  ohe  only  ionized  constituents  and  that  no  neg«» 

ative  ions  are  present,  so  that  the  atmosphere  consists  of 

+ 

neutral  particles  and  a  neutral  Anlxture  of  0  ,  He  ,  and  H  , 
and  electrons  in  diffusive  equlllbriuni.  The  ions  will  be  re¬ 
ferred  to  by  subscripts  1  and  electrons  by  subscripts  e. 

It  is  assumed  that  the  iines  of  force  act  as  barriers  across 
which  charged  particles  cannot  flow.  Iz  is  instructive  to  con 
sider  the  variation  of  partial  pressure  along  a  tube  of  force. 
We  can  then  write 


dp^  =  -m^n^f  ds  -  n^eEds 
e  e  e  e 

dp^  =  -m^n^f  ds  +  n^eEds 


1  (3, 


where  the  llsld  F.  (Fig.  3)  arises  from  she  small  separation 
between  charges  due  to  the  fact  that  electrons  are  lighter 
than  the  positive  ions  and,  therefore,  tend  to  move  upwards 
iKsnge  (i960)]. 


p,+<Sp, 

pi*<«Pi 


FIG.  S.  TO  ZU4JSTIUTE  AN  EUaCETT  OF  VOLUME  TAKEN  AL(H<0  A  LINE  OF  FORCE. 
Tlte  pMrtlitl  pr«ssur«8  of  th«  ions  and  olectrons,  p^  and  p  ,  and  the 
electric  field  E  (Mange  IS8O)  are  also  shown.  * 


SEL-63-110 


-  8  - 


2 


In  these  equations,  m^  and  m^  are  the  mass  of  the 
electron  and  of  the  1th  Ion,  n^  and  are  the  electron 

and  1th  ion  densities,  e  is  the  electronic  charge,  and 
and  represent  the  partial  pressures  of  the  electrons  and 

the  1th  ions,  respectively.  It  is  assumed  (Pig.  3)  that  the 
pressure  changes  over  a  small  volume  contained  within  the  tube 
of  force  by  an  amount  dp  over  the  distance  ds  measured  a- 
long  the  line  of  force.  The  direction  of  the  electric  field 
E  is  as  shown  in  Fig.  3.  It  should  be  noted  that  the  gradi¬ 
ent  of  partial  pressure  changes  direction  according  to  whether 
the  point  under  consideration  is  inside  or  outside  the  lines 
marked  CD  in  Fig.  2.  From  Dalton* s  law  for  the  partial  pres¬ 
sures  of  the  constituents  of  a  gas  in  diffusive  equilibrium, 
we  know  that  each  constltuf^nt,  in  this  case  each  ion  species, 
and  the  electrons,  has  a  partial  pressure  corresponding  to  that 
which  would  have  if  it  were  the  only  constituent  present. 


Ke  .T'  then  write 


^e  e  e 


Pi  “  "i'^i 


and  refer  to  the  electron  and  Ion  temperatures  re¬ 

spectively,  and  in  the  subsequent  theory  it  is  assumed  that 
thesr  can  be  different.  Since  the  atmosphere  at  every  point 
Is  ta.  ;n  to  be  electrically  neutral,  we  may  write 


n  =  r  ni 
®  i  ^ 


which  leads  to 


so  that 


E  jfr* 
1  ^1 


dp, 

1^ 
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Prom  Eqs.  (5)  and  (9)  the  electric  field  can  be  calculated: 


-  fds  -  eE  ds  =  t  -  fds  +  eE  T  6 


Z.  KiO 

i  -^1 


®^^1 


eE  = 


^  -n 
1  ^1 


""e  ^  ^  ’"1 

ST"  +  ^ 

•^e  1  -^1 


.  f 


eE  = 


I  iq  “l"!  -  ”e"e 

- '^e - 

"e  ^  "1 


.  f 


Since  ^  I850  and  ^  1  the  last  expression  can 

be  simplified  to  give 


Z  JR—  n.  ra. 

,  Ti  1  1  ra 

eE  =  i — i-w -  .  f  =  ^  f 

r,  4-^  e  „ 

*^e  1 


(10) 


where  ra^  is  the  "temperature-weighted  ion  mass  average" 


m,  =.  2 


Z  jfT’  n,m. 

1  ^1  ^  ^ 

n  't*  Z  m  Oj 
''  11 


(11) 


Using  (10)  In  (5)1 


m. 


<iPe=  f 


/  ®+ 

dPi=  -nj^Cm^-  ^)  f  ds 


(12) 
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From  (ll).  If  TgAi  ^ 


5Z  ^ 

2  ^  ^  m .  ^  y  p 

i  ^  i  ^ 


and  thus  we  may  write  (12)  as 


‘^Pe  =  ■  "e  2"  ^ 


“+ 

dp^  --  -  “  ^)  ^ 


(13) 


From  Eq.  (6)  we  may  hen  write 


dPg  =  kd(ngTg) 


dpj^  =  kd(nj^Tj^) 


(I't) 


d(n^T^) 


m, 

cl(niTi)  (mj^  -  ^) 


(15) 


Equations  (15)  have  a  number  of  Important  consequences,  which 
are  discussed  below.  It  should  be  noted  that  If  a  constant 
temperature  Is  assumed  along  a  line  of  force  for  both  the 
electrons  and  ions.  Equations  (15)  show  that: 

1.  Since  m.  Is  always  greater  than  zero,  the  electron 
density  decreases "With  s  and  thus  also  with  height. 

2.  If  m.  >  ra./2,  the  Ion  density  n.  decreases  with  s 
and  thus  with  height.  This  condition  is  always  satis¬ 
fied  by  the  heaviest  constituent  and  also  by  a  light 
constituent  if  it,  or  a  still  lighter  one  is  strongly 
predominant  (see  Fig.  4,  page  21). 
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3.  If  nu  <  m./2,  the  ion  density  n*  Increases  with  s 
and  thus  with  height.  This  happens  with  a  light  constlt 
uent  In  the  presence  of  a  predominant  heavier  one. 

4.  If  nu  »  m./?,  the  Ion  density  of  constituent  1  (can¬ 
not  he  the  heaviest)  passes  through  a  maximum. 

These  equations  are  basic  and  contain  the  essential  in¬ 
formation  about  the  distribution  of  both  electron  and  the  Ion 
densities  with  height.  In  order  to  solve  them  to  obtain  this 
distribution,  however.  It  Is  necessary  to  make  a  number  of 
changes  of  variable  to  simplify  the  mathematics.  For  conven¬ 
ience,  distances  along  a  line  of  force  will  be  measured  from 
a  reference  level  at  a  distance  r^  (Pig.  l)  from  the  center 
of  the  earth.  All  quantities  measured  at  this  level  will  be 
indicated  by  the  subscript  zero.  Integrating  Eqs.  (15)  be¬ 
tween  the  limits  s  =  0  and  s  =  s‘,  the  electron  and  Ion 
densities  at  the  distance  s’  from  the  reference  level  along 
the  line  of  force  are  given  by 


(16) 


In  which  the  ion  and  electron  temperatures  must  be  given  as 
functions  of  s . 

In  the  above  Integrals  it  Is  convenient  to  separate  the 
effects  of  T  and  f,  as  follows; 


r®’  ”+80  ff  %o] 

J„  sErrlsT  T— 

o  eo  1^0  e 


ds 


ni.dz 

+ 


e 


(17) 
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III.  m, 

|.s'(ra^  -  -  — )So  ff 

Jo  '  I  K  ^j 


lo  o 


f  ^  "*+ 

J  (m^  -  5-)dEj^ 


(17) 

(Cont'd) 


In  these  equations,  and  z^  are  the  "temperature- 

modlfied  geopotential  heights"  for  electrons  and  each  ion 
species  respectively  and  are  defined  by 


T 

_  f  -"eo 
'  So 


T 

d'^  -L.tl2.6s 

”  So 


(18) 


from  which  we  get,  setting  =  0  at  the  reference  lev¬ 

el  s  =  0, 


~  j  cf  T  ( s ) 

^s=0  ^e'  ' 


'1  =  Js=o5;^'- 


(19) 


Thus,  if  Tgo'^e^®^  different  from 

then,  at  a  point  s’,  the  quantities  z^  and  z^  have  differ¬ 
ent  values. 

Equations  (l6)  can  now  be  written  in  terras  of  the  temper¬ 
ature-modified  geopotential  heights  z^  and  z^^  as: 


=  "eo  2l?r^ 
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“  "10  f^Ts^  ®*pI-  B^l  •  ®*p(5i?rf;  V®il 


(20) 

(Cont’d) 


Where 


H 

"l  -  miSo 


(20a) 


is  the  scale  height  of  the  non-ionlzed  atomic  species  1,  at 
the  reference  level. 

In  Eq.  (20),  n^(s‘)  and  n^^Cs*)  are  given  In  terms  of 

Integrals  In  [given  by  Eq.  (ll),  and  still  unknown]  with 

respect  to  and  z^.  To  eliminate  a  relationship  must 

be  established  between  and  z^.  Equations  (19)  show  that 
this  can  be  accomplished  either  by  knowing  and 

as  functions  of  s  or  by  assuming 


eo  lo  4.  f 


(21) 


This  means  that  the  electrons  and  all  the  ions  have  the  same 
ten5>erature  variation  with  s,  the  distance  along  a  line  of 
force.  If  Eq.  (2l)  applies,  then  we  may  write  (19)  as 


s«s 


=  z^  =  z  =  J  t(s)  ds 


s=o  ®o 


(22) 


Thus,  the  electron  and  ion  densities  can  be  written  In 
terms  of  the  same  z  as: 


eo  o 


(23) 
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ni(s')  =  niot(s')  exp  (-  |-]  expf^^  m^dxj 


(23) 


Because  of  charge  neutrality. 


n^(s')  =  t  =  ^plot(  =  ’)  ">+H  1(2“) 


The  last  exr'Jnentlal  can  be  calculate(3  from  the  first  of 
Eqs.  (23): 


■io"o 


Sc 

io  ^eo"o 


=  H2h;;Cv^1)  =(^ 


®°  t(s') 


1“ 


(25) 


where 


T 

n  eo 

Ci  =  ^ — 

^  ^io 


Equation  (24)  becomes: 


n_(s*)  =  z\ 
®  1 


n^gt(s')  exp(-  I-)  j_^t(s-)j 

1  Q 


or: 


n  (s») 

-2 -  =  -L-  z 

%o  ^eo  1 


l+C, 


n  ^1 
eo 


n,„t(s.)  ^  axp(-  |-)(5^) 


(26) 


(27) 


Normalizing  the  ion  densities  at  the  reference  level  to 
the  value  for  the  heaviest  ion  (O^),  we  write 


^  ^  =  Pi  =  ^1^^2^^3 
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wherv 


1^  =  1 


^2  ~  ^2o/^lo 
^3  "  ^3o/’^lo 


(28) 


so  that  T]  is  the  electron  density  divided  by  the  oxygen  ion 
density  at  the  reference  level. 

Equation  (27)  becomes 


IV  (s’)  ,  I  1+C.  _  n„  C. 

"  exp(-  1^)  (5^)  ^ 


(29) 


in  this  equation,  the  three  first  factors  in  the  summa¬ 
tion  are  known,  and  the  equation  must  be  solved  for  (ng(s’)/ 

h«^) . 
eo ' 

Suppose  now  that  the  ratios  of  ion  temperature  to  elec¬ 
tron  temperature  at  the  reference  level  are  the  same  for  all 
ions,  that  is: 


G  —  —  G 

i  -  Tlo 


1  =  1,  2,  3 


(30) 


Then  Eq.  (29)  becomes 

n^(s') 

(-^— ) 

^eo 


=  t(s') 


Ti  f  l^^i  ®^P("  j 


from  which  we  get,  finally,  the  electron  density 

'^eo  (1  _  f.  z  vl  1 

=  57^7  If  [  «P(  -H^)|  1 


(31) 
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From  the  second  of  Eqs.  (23)  and  (25)  we  get  the  Ion 


densities ; 


n.(s')  T  C+1 

^  exp( 


(32) 


The  ion  density  may  also  be  referred  to  the  electron  density 
at  the  reference  level  by  Eqs.  (28)  so  that 


n.  (s')  n.  T  ^  r  n 

n-.  '  “  57*  [nTs°)’] 


(33) 


Or,  using  Eq.  (31 ), 


“  "n  T  (s' )  " 

6 


|-)  [s - ^  ^  (3^1) 

1  1  exp(-  jq-) 


The  distribution  of  electrons  and  Ions  in  the  exosphere  are 
given  by  the  Eqs.  (31 )  and  (33)  in  terms  of  the  variable  z  , 
the  temperature-modified  geopotential  height  defined  by  Eq. 
(22).  These  equations  provide  us  with  the  information  we  re¬ 
quire  to  give  the  distribution  of  electrons  and  ions  along  a 
line  of  force,  and  hence,  vertical  profiles.  It  is  shown  In 
Section  V  that  these  general  solutions  reduce  to  the  equations 
given  by  Bauer  (1962)  for  the  particular  cases  he  considers. 
The  derived  distribution  Eqs.  (3l),  (33),  and  (3^)  are  dis¬ 
cussed  in  detail  below. 


C.  FORM  OF  SOLUTION  FOR  EQUILIBRIUM  ION  AND  ELECTRON  DENSITIES 

For  the  case  in  which  the  ion  and  electron  temperatures  are 
equal  at  the  reference  level  [and  hence,  elsewhere,  by  Eq.  (21)], 
we  have  C=1  and  Eq.  (31 )  becomes 


(35) 
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If  the  electron  temperature  Is,  say,  twice  the  ion  temperature. 


OillU 


ne(s’) 


n 


eo 


'  ijjl 


1 

T 


It  Is  clear  that  If  such  a  difference  between  Ion  and  electron 
temperatures  occurs,  then  the  electron  densities  may  be  changed 
by  a  factor  of  10  at  very  high  altitudes. 

If  the  temperature  along  a  line  of  force  is  constant  so 
that  TgoAg(3’)=l^  then 


(36) 


Prom  Eq,  (34)  we  have  the  ratio  of  densities  of  two  ions: 


n. 


-  H, 


) 


T1  H,  -  H, 

^exp(z  -fqnp) 


(37) 


The  Ion  densities  will  be  equal  at  a  value  of  z  =z 


defined  by 


z  in  h 

Hi  =  1"  < 


IJ 


(38) 


Since  the  subscripts  1,  2,  and  3  refer,  respectively,  to 
+  +  + 

0  ,  He  ,  and  H  ,  the  scale  heights  are 


Hg  -  4 

=  16 


(39) 
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Hence,  the  values  of  z  at  whloii  n-=n^, 

n2=n^  are,  respectively,  from  equations  (38): 

He*^)  =  -  ^  >  0  if  \  <  1 

z^3(ot-^  H'^)  =  -  ^  H^lnTI^,  >  0  if  11^  <  i 


Z23(Het^H'^)  =  "  3  H^ln  >  0  if  <  T,^ 


(40) 


Thus,  two  Ions  will  have  the  same  density  at  some  height 
above  the  reference  level  if,  at  the  reference  level,  the 
heavier  one  is  more  abundant.  Also,  from  Eq.  (37),  at  z 
larger  tnan  ,  the  ion  i  will  have  higher  density  than 
the  ion  J  if  it  is  lighter.  Due  to  this  fact,  the 
defined  by  Eqs.  (40)  are  called  transition  values  of  z, 
and  these  are  indicated  in  the  same  equations. 

It  is  interesting  to  consider  the  conditions  for  which 

4-  4-  +  + 

the  transition  He-^H  happens  above  the  transition  0-^He  . 
From  Eqs.  (40)  this  condition  is 


Ti^  >  .iX 


(iti) 


The  distribution  of  electron  density  as  a  function  of 
(31)  may  thus  be  written  [using  (39)] 


ng(s') 


eo  rexp(-x  )+T\  exp(-x  /4)+ll-.exp( -x  /I6)  il+O 

JTJ  1 - ii - ^ - 1  (^2) 


where 


X  = 


(43) 
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1. 


It*  the  reference  level  is  low  enough  so  that  the  0^ 
is  strongly  predominant  there,  then 


^2  1 
«  1 


I  (44) 


For  low  values  of  x  the  exponentials  in  (42)  are  close 
to  unity,  and  we  have 


n,(s')  ^  T 


eo 


exp(-  ^^) 


(45) 


Thus,  apart  from  the  factor  '^eoA«(s’)#  the  electron 
density  varies  exponentially  in°  z,  wlt:h  scale  height 
(1+C)H^. 

2.  For  large  values  of  x  the  two  first  exponentials  in 
(43)  become  much  smaller  than  the  third,  so  thac  they 
can  be  disregarded  and 


Apart  from  the  factor  the  electron  density 

presents,  as  a  function  Sr  f,  an  exponential  behavior, 
with  scale  height  (l+C)  H^,  and  an  asymptotic  value 

1/(1+C) 

(^) 

for  z=G , 

For  the  case  C=l,  the  electron  and  ton  distributions 
calculated  as  a  function  of  z  from  (3l)  and  (33)  respec¬ 

tively  are  shown  in  Pig.  4,  for  T=1000^K,  112^0.02  and 
'T’2-0.0016  [Bauer  (1963)].  The  dependence  of  n^  and  n^ 
on  altitude  is  implicit  in  these  curves  of  Fig.  4  through  the 
dependence  on  z  (see  Pig.  5#  p.  25).  A  discussion  of  the 
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z,  km 


3500- 


RELATIVE  DENSITIES 


FIG.  4.  DISTRIBUTION  OF  ELECTRCW  AND  I(H<  DENSITIES  WITO  THE 
TEMPERATURE-MODIFIED  GEOPOTENTIAL  HEIGHT,  z,  FOR  THE  CASE 
f'  =  1  AND  OONSTANT  TEMPERATURE  T  =  1000  ‘’K. 
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resulting  electron  and  ion  density  distributions  as  a  function 
of  height  is  deferred  to  Section  IV.  The  curves  of  Pig.  4, 
however,  are  quite  fundamental  to  an  understanding  of  the  pro¬ 
cesses  Involved  in  the  theory  of  the  exosphere  and  the  curve 

for  n  as  a  function  of  z  is  the  same  as  that  given  by 
e 

Bauer  (1963).  The  quantity  z  is,  however,  in  the  present 
work  given  a  different  meaning  from  that  used  by  Bauer,  since 
it  contains  a  dependence  on  temperature  and  centrifugal  force. 
The  resulting  differences  predicted  by  the  theory  presented 
here  compared  with  the  results  given  by  Bauer  (1963)  and  others 
are  discussed  in  Section  V, 

Using  Eq.  (33)  it  is  seen  that,  for  C=l, 


Hence,  the  geometrical  mean  of 

quantity  In  parentheses.  For  an  isothermal  exosphere,  this 
quantity  happen ««  to  have  a  linear  variation  with  z  ,  when 
plotted  on  semi-log  paper,  with  "scale  height”  and  pass¬ 
es  through  the  point  at  z  =0. 

This  important  property  was  used  to  obtain  the  ion  den¬ 
sities  from  the  electron  density,  in  Pig.  4.  More  details  are 
given  in  Appendix  B, 

The  main  source  of  difficulty  in  obtaining  the  solutions 
of  (31)  and  {34)  is  to  be  found  in  the  derivation  of  the 
quantity  z  .  This  is  discussed  in  the  next  section. 
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D.  CALCULATION  OF  TEMPERATURE-MODIFIED  GEO  POTENTIAL  HEIGHT  2 
It  will  be  recalled  that  the  quantity  z  is  defined  by 


,.s'f(e  ,s)  T 

^  (®o-s')  =  J.  — r — TT 


CO 


('*7) 


In  order  to  evaluate  z  it  is  necessary  to  transform  the 
variable  of  integration  in  the  above  equation  from  s  to  9 
where  6  is  the  geomagnetic  latitude  of  the  point  concerned 
on  the  line  of  force.  In  Appendix  A  it  is  shown  that  the 
element  of  distance  ds  along  the  line  of  force  may  be  writ¬ 
ten  in  terms  of  the  subtended  angular  increment  d9  by  the 
relationship 


as  =  — 2- —  cos  e  vi+3  sln^  0  d0 

cos^ ®o 


(48) 


Substituting  for  ds  from  (48)  and  for  f  from  (4)  in  the 
above  expression  for  z  and  using  (a4).  Appendix  A,  we  have 

z  =  -J  s^flr  — cos  e  /l+3sln^9  d9 

''9  -^e'  '  008^^9 

o 


.9_  T 


2 

cos  9. 


=  1  rUT  ”0  ^  cose  - 

'^9*  ^e'®^  °  *  cos^Q 


^  ^o  cos^9 


^o  cos  9. 


COSY  Id 9 


(49) 


where  cos0  is  given  by  (3a)  and  cosy  by  (3b)  (Fig.  l) 


Thus 


,2cos^9^  3n^r„cos^9, 

=  I  ^  ^o  rf%  V  ■  ■  - ir - % -  Sin0d9 

•>6  °  V  n-4tan  9  ‘  cos  9  g^cos  9  * 

®o  o 


r2eos^9^  3n^r^cos59  , 

'  J  rffr  — 3-^ - 1 

fit  ^  cr  r*ri«  Q 


cos-^S 


g^cos  9^ 
o  o 


(50) 


-  23  - 


SEL-63-110 


since 


9  ^  , 

l+4tan  6 


This  expression  for  z  ,  the  temperature -modified  geopoten- 

tlal  height,  If  substituted  In  Eqs.  (3l)  and  (3^)  will  give 

the  density  distribution  of  electrons  and  ions  along  a  line 

of  force  for  any  given  temperature  distribution  along  a  line 

of  force.  In  general,  the  equation  will  have  to  be  computed 

numerically  on  a  digital  computer.  It  is  Interesting  to  note, 

however,  that  In  the  case  of  a  constant  temperature  along  a 

line  of  force  so  that  T  (6)  =  we  have 

e  eo 


z  =  2r'rt®os^e  X  .  d9  - 


COS-^Q 

e. 


O  2  ^2 


° 


=  2rQCos 


°^2oos®9^e' 


eo®°®  ®o 


toos^el  ° 


sln0cos'^0d9 


9. 


=  r. 


cos^9 

I 

cos  9 


2 

n^r 


26, 


6  ’ 

o  f 2«  cos  6 

_  ^cos  e - If— 

^  COS  9_ 


(51) 


SO  that  In  the  simple  case  of  an  isothermal  atmosphere  i:. 
has  an  analytical  solution. 

The  resultant  variation  of  z  with  altitude  is  shown 
In  Fig.  5  for  different  latitudes.  If  the  centrifugal  force 
Is  disregarded,  z  is  a  function  of  altitude  only  and  is 
Independent  of  latitude  (curve  at  right).  It  approaches  a 
limiting  value  (equal  to  the  distance  from  the  reference  level 
-to  the  center  of  the  earth)  as  the  altitude  goes  to  infinity. 
Inclusion  of  the  centrifiigal  force  implies  a  decrease  in  z  , 
especially  at  high  altitudes  and  low  latitudes.  The  quantity 
z  then  has  a  maximum  value  at  approximately  6  earth  radii 
for  low  latitudes  and  approximately  8  earth  radii  for  higher 
latitudes. 
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GEOPOTENTIAL  HEIGHT,  Z,  km 


5.  VARIATION  OF  TEMPERATURE-MODIFIED  GEOPOTENTIAL  HEIGHT, 
WITH  ALTITUDE,  h.  F(Ht  AN  ISOTHERMAL  EXOSPHERE. 
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It  should  be  noted,  of  course,  that  In  Eqs.  (31)  and  (34) 
it  is  possible  for  the  temperature  along  different  lines  of 
force  to  be  different  and  to  be  distributed  along  each  line 
of  force  in  a  different  way.  As  stated  above,  the  Eqs.  (3l) 
and  (34)  in  conjunction  with  the  general  relationship  (50)  may 
be  used  to  compute  the  profile  of  electron  density  along  a  line 
of  force.  If,  however,  we  require  only  the  ratio  of  electron 
densities  in  the  northern  and  southern  hemispheres  at  the  base 
level  (see  section  E  below),  then  it  is  useful  to  consider  the 
form  of  Eq.  (51)  when  9*  is  small  so  that  the  point  (9^,9*) 
lies  near  the  equatorial  plane.  For  9»-0  the  value  z  of 
z  from  (51)  becomes 


o 

(52) 

E,  CALCULATION  OP  RATIO  OP  ELECTRON  DENSITIES  AT  CONJUGATE 
POINTS 

Attention  will  be  given  now  to  a  rather  interesting  con¬ 
sequence  of  the  theory  presented,  namely  the  ratio  R  that 

eo 

must  exist  between  the  electron  densities  at  the  reference 
level  at  two  magnetically  conjugate  points  to  satisfy  conti¬ 
nuity  at  the  top  of  the  field  line. 

It  will  be  assumed  that  the  ions  and  electrons  have  the 
same  temperature,  and  hence  C=1  in  Eq.  (30). 

1.  Temperature  Constant  in  Each  Hemisphei^e 

As  a  first  approximation,  suppose  that  the  average 
temperature  in  the  southern  (Summer)  hemisphere,  Tg,  is 
greater  than  the  average  temperature  in  the  northern  (Winter) 
hemisphere,  Tj^,  (Pig.  6)  so  that 


[l-eos^So]  + 


2 

COS  9^- 

o 


4 

cos  9 


with  K  a  constant  greater  than  1, 
SEL-63-110  -  26  - 


(53) 


FIG.  6.  PARAMETERS  USED  TO  CALCULATE  R«o  AT  TEE  REFBHEMCS  LEVEL  AT 
TWO  MAGRETICALLY  CONJUGATE  POINTS  VBEN  THE  TEMPERATURE  IN  THE  MORTHERN 
AND  SOUTHERN  HEMISPHERES  ARE  DIFFERENT  BUT  COMSTANT  IN  EACH  HBUSPBERE. 


Given  the  electron  density  n^^  (and  Ion  densities, 
through  ^2  ^3)  ^  reference  point  on  a  field  line.  It 

is  possible  to  calculate  the  electron  density  at  any  other 
point  (defined  by  0^,9')  on  the  same  field  line,  using  Eqs. 

(36)  and  (51). 

Let  A  be  the  ratio  given  by  Eq.  (36)j 

it  is  clearly  a  function  of  temperature,  through  the  values 
of  and  Thus,  for  the  northern  and  southern  hemi¬ 

spheres  respectively. 


where  is  the  value  of  z  for  e»-  0  and  Is  given  by 

(52). 
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As  the  electron  pressure  must  be  the  same  In  points 
1  and  2  of  Pig.  6: 

"eN'^N  =  Pn  =  Ps  =  "es'^S  <55) 


Prom  (53),  (5^)  and  (55), 


n 


R. 


eoN 


(ej 


eo 


n 


eoS 


ro 


KjT 


(56) 


This  expression  permits  the  calculation  of  the  elec¬ 
tron  density  at  the  reference  level  in  the  northern  hemisphere 
if  the  electron  density  at  the  same  level,  at  the  conjugate 
point,  is  known.  It  is  useful  to  make  the  following  approxi¬ 
mations  , 

Prom  Pig.  4  it  is  seen  that  for  z  greater  than  a- 
bout  2000  km  (corresponding  to  at  6^  >  35°),  Eq.  (46) 
holds  so  that  with  the  present  assumptions,  using  Eq.  (20a), 

Ajj  and  Ag  are  given  by 


Ag  w/n^g/Tig  exp(-Zj^/2H^g) 


(57) 


Thus,  at  latitudes  greater  than  about  35°,  Eq.  (56)  has  the 
approximate  form 
n. 


R 


eo  n 


*eoN^_®o^ 
eoS^®o^ 


4. 


exp(- 

3NS 


W.. 


n-) 


3N 


(58) 
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The  corresponding  equation  given  by  Rothwell  (1962) 
may  be  derived  as  a  particular  case  of  the  more  general  ex¬ 
pression  (58)  in  which  the  dependence  of  m,  g  and  T  on 
altitude  may  be  properly  allowed  for. 


In  the  preceding  section,  it  was  assumed  as  a  conven¬ 
ient  approximation  that  there  was  a  discontinuity  of  tempera¬ 
ture  between  two  points  Just  north  and  south  of  the  equatorial 
plane  (Fig.  6).  For  a  complete  treatment  it  is  better  to  pro¬ 
ceed  as  follows. 


FIG.  7.  PARAMETERS  USED  TO  CALCULATE  R«o  AT  THE  REPERBRCB  JJSVEL  AT  THD 
MAGWETICALLT  CONJUGATE  POINTS  WHEN  THE  TEMPERATURE  IS  TAKEN  AS  A  CON¬ 
TINUOUS  FUNCTION  T]j(8)  in  THE  N(XtlHERN  HBU SPHERE,  AMD  T3(s)  IN  1BI 
SOUTHERN  HEMISPHERE.  At  points  1  and  2,  tba  toaparaturas  ara  tha  aaaa. 


Figure  7  indicates  the  temperature  as  a  function  of 
distance  along  a  line  of  force.  The  temperature  must  change 
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along  the  field  line  since  it  Is  supposed  that,  at  the  refer¬ 
ence  level,  it  has  different  values  In  the  northern  and  south¬ 
ern  hemispheres 


(59) 


with  K  a  constant  greater  than  1. 

The  distance  s  is  clearly  a  function  of  8^  and 
0,  and  we  call  the  value  of  s  corresponding  to 

0=0  (details  are  given  in  Appendix  A). 

The  functions  Tjj(s)  and  Tg(s)  are  such  that  the 
temperature  is  a  continuous  function  of  s,  including  the 
point  Sj^,  that  is: 


The  assumed  form  for  'tjj(s)  and  Tg(s)  is  discussed 
in  Section  IV. 

As  the  temperature  is  now  a  function  of  s,  Eq.  (50) 
must  be  used  to  calculate  z  ,  and  different  values  will  re¬ 


3 

and  4, 

Fig.  7: 

T 

^oN 

r2cos^9_ 

3Q^r 

0  cos-'01 

Je- 

^  cos-^e 

g  4  J* 

®0  cos  8 

de 


-  '•0J3,  '■L^Z3 


Sin  ef 


2cos'^e^  3n^r  5 
0-^0  co'’-^ 


COS”^0 


cos 


(61) 


These  equations,  together  with  Eq.  (35),  give  the 
electron  density  distributions  in  the  northern  and  southern 
hemisphere o. 


n 


eoN 


■'oN 

rrnr 


V  (62) 
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nes(s') 


=  Tirin'  1  [’’is  jqi’l  1 


(6a) 


For  s*=  the  electron  densities  are  the  same  in 
the  northern  and  southern  hemispheres  (points  1  and  2,  Pig.  ?)> 
and  so  are  the  temperatures.  Thus,  the  last  equations,  togeth¬ 
er  with  (60)  give 


i  h’lS  }  "F 


'S  ^  -IN 


(63) 


where  and  Zj^g  are  given  by  the  corresponding  Eqs.  (6l) 

with  6'  -  0.  For  sufficiently  large  values  of  and 

1*10 

z^,  such  that  Hq.  (46)  holds,  Sq.  (63)  may  be  written  as 


1 

-  J  N  3ST 

p  ~  exp 


(^4) 


in  which  use  has  been  made  of  Eq.  (20a)  and  (59). 
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IV.  RESULTS  OP  THE  THEORY 


A.  INTRODUCTION 

In  this  section,  relative  electron  and  ion  density 
alstributlrns  calculated  on  the  basis  of  the  theory  outlined 
in  the  previous  sections  are  presented.  The  results  are 
given  as  distributions  calculated  as  functions  of  the  tempera¬ 
ture  modified  geopotential  height  z,  as  functions  of  the 
distance  s  measured  along  a  line  of  force  and  also  as  functions 
of  the  altitude,  h.  These  distributions  are  referred  to  as 
N(2),  N{s)  and  N(h)  curves  respectively.  The  computations 
are  based  on  a  number  of  different  assumptions  about  the 
exospheric  temperature  and  about  the  ionic  composition  of  the 
atmosphere  at  500  kilometers.  In  particular,  the  distributions 
are  derived 

a)  when  a  constant  temperature  is  assumed  along  a  line 
of  force 

B)  when  the  temperature  is  assumed  to  vary  with  distance, 
s,  along  3  lire  of  force. 

In  addition,  two  sets  of  assumptions  are  made  about  the 
composition  at  the  base  level  usually  taken  to  be  500  kllometen 
above  the  earth’s  su^'face.  These  are 

1.  that  the  composition  at  500  kilometers  is  independent 
of  the  tei^perature  at  that  level,  and 

2.  that  the  composition  at  500  kilometers  is  a  strong 
function  of  the  temperature  at  that  level. 

The  N(h)  calculations  referred  to  above  were  carried  out 
both  when  a  proper  allowance  was  made  for  the  effect  of  the 
centrifugal  force  arising  from  the  earth’s  rotation  and  also 
when  this  force  was  neglected.  Subsequently,  use  is  made  of 
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observational  data  from  the  Alouette  satellite  (Thomas  and 
Sader  1963)  to  convert  the  relative  electron  density  profiles 
into  absolute  profiles  based  on  the  electron  density  observed 
at  approximately  1000  kilometers.  This  leads  to  a  set  of 
theoretical  profiles  calculated  using  reasonable  experimentally 
observed  magnitudes  for  ion  composition  (Bauer  1963),  tempera¬ 
ture  (Harris  and  Prelster  1962),  and  electron  density  at  lOOG; 
km  (Alouette  observations ) .  These  theoretical  profiles  are 
then  compared  vfith  a  number  of  whistler  and  other  observations 
of  electron  density  in  the  earth’s  exosphere.  Tables  i  and  2,pg. 
61, 62, give  a  summary  of  the  way  in  which  the  calculations 
were  developed  and  in  which  the  results  are  presented.  Further 
details  are  noted  in  tne  ensuing  sections. 

B.  THE  ELECTRON  AND  ION  DISTRIBUTIONS 

The  electron  and  ion  distributions  presented  below  were 
calculated  using  a  digital  computer  program.  The  N(2)  curves 
were  computed  from  equation  (35)  using  the  assumption  that 
the  temperature  was  constant,  the  N(s)  curves  were  computed 
from  equations  (35),  (51)  and  (A4),  the  N(h)  profiles 
being  subsequently  computed  from  the  N(5)  profiles  using 
electron  density  observations  at  1000  kilometers  as  measured 
by  the  Alouette  satellite. 

The  calculations  in  which  the  temperature  is  assumed 
constant  along  the  line  of  force  of  the  earth’s  magnetic 
field  are  denoted  by  the  letter  A  (Table  1,  pg.6l).  For  each 
value  of  the  temperature  assumed,  it  is  necessary  to  specify 
the  composition  at  the  base  level  at  500  kilometers.  This 
was  taken  from  the  data  presented  by  Bauer  (196.3)  foJ’  the 
helium  and  hydrogen  ion  densities  relative  to  the  oxygen 
ion  density  as  a  function  of  the  temperature.  The  relevant 
equations  are 
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TIgCTo)  =  0.2  X  10 

^3(Tq)  =  0.16  X  10 


(65) 


In  this  equation,  is  the  temperature  assumed  at  the 
reference  level.  The  variation  of  Tl^  and  with  temperature, 
as  given  by  equations  (65}  is  shown  in  Pig.  8. 

The  exospheric  temperature  enters  the  calculations  in 
two  ways.  Firstly,  via  the  temperature  modified  geopotential 
height  z,  and  secondly,  via  its  effect  on  the  composition 
assumed  at  the  base  level  through  equations  (65).  It  is, 
however,  interesting  to  carry  out  computations  in  which 
the  composition  at  500  kilometers  is  held  constant  neverthe¬ 
less  allowing  the  temperature  to  vary.  In  this  case  the 
values  for  and  Tj^  were  taken  from  equations  {65)  using 
T^  =  1000®K  and  these  results  are  referred  to  as  "Composition 
1,"  Table  1. 

The  calculations  in  which  equations  (65)  were  applied 
to  give  the  relative  densities  at  500  kilometers  are  referred 
to  as  "Compositions  2a,b,c,"  etc,  depending  on  the  value  of 
the  temperature  at  500  kilometers  (see  Table  1).  In  other 
words,  in  "Compositions  2,"  once  the  temperature  at  the 
base  level  is  specified,  the  relative  ion  densities  are 
automatically  specified  at  that  level  by  equations  (65). 

Some  calculations  are  also  performed  using  the  assumption 
that  the  temperature  is  changing  along  the  field  line. 

These  results  are  denoted  by  the  letter  B  (Table  1)  and  will 
be  described  later. 

Details  relating  to  the  curves  shown  in  Pigs.  9  -  18 
are  summarized  in  Tables  1  and  2  which  precede  the  figures 

(pg.  61  ). 
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In  Figures  9-1  -  e  (see  Tables  1  and  2),  electron  and 
ion  densities  are  plotted  as  functions  of  the  temperature- 
modified  geopotential  height  z  for  different  values  of 
the  (constant)  temperature  (assumed  independent  of  s )  and 
for  different  compositions.  The  distributions  shown  in 
Figs.  9a  -  e>  ho  not  show  marked  "kinks"  or  cusps  and  the 
Ion  transition  levels  do  not  show  up  clearly  in  the  electron 
density  curves.  These  curves  are  discussed  in  greater 
detail  in  Appendix  B. 

The  electron  densities  of  Figs.  9a  -  e  are  summarized 
in  Pig.  10,  where  the  relative  effects  of  temperature  and 
composition  changes  can  be  seen.  At  low  values  of  z,  where 
0^  is  predominant  for  all  curves,  the  slopes  are  proportional 
to  temperature.  As  we  go  to  higher  values  of  z  the  effect  of 
the  Ion  composition  at  the  base  becomes  relatively  more 
Important  and  the  slopes  are  no  longer  simply  proportional 
to  temperature.  For  Instance,  at  z  =  3000  km,  the  slope 
of  curve  e  (2000°K)  is  less  than  the  slope  of  curve  b 
(1000°K).  This  is  sc  because,  although  for  curve  b  at 
z  =  3000  km,  the  is  strongly  predominant  (Fig.  9d)*,  this 
is  not  the  case  for  curve  e  (Pig.  9e). 

It  is  seen  from  the  curves  of  Pigs.  9a  -  e  that  the 
ion  transition  levels  (the  levels  at  which  n(0'*’)  =  n(He'*^), 
etc.),  change  with  temperature.  The  variation  of  the  ion 
transition  levels  with  temperature  is  shown  in  Pig.  11  for 
two  different  assumptions  about  the  composition  at  the  base 
level,  namely,  that  the  composition  is  independent  of 
temperature  (continuous  lines)  and  that  the  composition  is 
temperature -dependent  as  indicated  by  equations  (65) 

(broken  lines).  An  Important  conclusion  which  can  be  drawn 
from  an  examination  of  Pigs.  93  -  e,  and  11  is  that  the 
transition  levels  approach  nearer  the  earth  at  night.  These 
curves  agree  with  those  given  by  Bauer  (I963). 
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The  N(z)  curves  may  be  transformed  using  equation  (51) 
into  curves  giving  N  as  a  function  of  0  and  9^  where  0  is 
the  angle  between  the  equatorial  plane  and  a  point  at  a 
distance  a  measured  along  a  field  line.  The  field  line  is 
assumed  to  Intersect  the  base  level  (500  km  above  the  earth) 
at  a  point  making  an  angle  0^  with  the  equatorial  plane 
(Pig.  1).  The  use  of  equation  (A4)  (Appendix  A)  converts 
the  N(0,0q)  curves  to  distributions  of  electron  density 
along  a  line  of  force  (N(s,0q)  curves).  Curves  of  this 
kind  are  given  in  Figs.  12a, b,  c>  for  9^=  45°,  55°  and  65° 
respei tively.  In  each  diagram  the  results  corresponding  to 
the  case  in  which  proper  allowance  has  been  made  for  centri¬ 
fugal  force  are  shown  as  continuous  lines.  The  corresponding 
results  for  the  cases  in  which  the  centrifugal  force  was 
neglected  are  shown  by  broken  lines.'  As  expected, 

the  effect  of  the  centrifugal  force  is  largest  at  the  higher 
latitudes,  and  tends  to  Increase  the  electron  density  at 
the  high  levels.  The  Increase  is  never  greater  than  about 
20^  over  the  range  of  heights  considered.  It  is  found  also 
that  along  the  field  line  at  45°  the  centrifugal  force  has 
no  appreciable  effect  (Fig.  12a).  The  effects  of  changes  in 
temperature  and  composition  are  seen  very  clearl"  in  Fig.  12a 
in  which  the  distance  scale  is  expanded. 

Similar  calculations  (denoted  by  the  letter  B,  Tables 
1  and  2)  were  carried  out  in  which  the  temperature  varied 
along  a  line  of  force.  The  temperature  distribution 
assumed  along  a  line  of  force  in  the  cases  labelled  B  is 
given  by  the  formula 

T(s)  =  T„  +  (T^  -  T„)  [:i%  (66) 

where  T^  is  the  temperature  at 't|ie  reference  level  (500 
kilometers),  T^  is  the  temperature  "St  the  top  of  the  field 
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xine,  D  Is  a  variable  distance  along  the  line  of  force 

4 

(taken  to  be  10  km),  and  the  distances  s  and  are 
defined  along  the  line  of  force  as  indicated  In  Pig.  7> 
page  29. 

In  order  to  be  able  to  calculate  the  ratios  of  electron 
densities  at  magnetically  conjugate  points  when  a  difference 
of  temperature  exists  between  hemispheres,  it  is  convenient 
to  apply  equation  (66)  separately  for  the  two  hemispheres. 
Thus  for  the  northern  (winter)  hemisphere,  and  were 
taken  to  be  1000° K  and  l800°K  respectively  and  the  resulting 
temperature  distribution  along  the  line  of  force  as  given 
by  (66)  for  the  northern  hemisphere  is  denoted  by  T^(s). 

For  the  southern  (summer)  hemisphere  =  1400°K  and 

Trj,  =  1800°K  and  equation  (66)  is  again  used  to  give  T2(s). 

The  distributions  Tj^(3)  and  ^2(3)  are  sketched  in  Fig.  13. 

Electron  densities  along  lines  of  force  were  calculated 
for  a  variable  temperature  along  the  line  of  force  using 
equations  (35)>  (50),  and  (a4).  The  results  are  shown  ip 
Figs.  l4a  -  d,  for  =  35®,  45°,  55®,  and  65®,  respectively. 

The  general  behaviour  of  the  curves  in  Figs.  14  is  similar 
to  that  of  the  curves  in  Pigs.  12.  For  instance,  an  increase 
in  temperature,  keeping  the  base  level  composition  constant, 
increases  the  electron  densities  in  both  cases,  and  a 
contrary  effect  is  produced  if  the  composition  is  supposed 
to  vary  as  given  by  equations  (65).  Also,  the  effect  of 
neglecting  the  centrii’ugal  force  is  negligible  below  8^  ?=  . 

It  is  im''ortant  to  note  that  the  calculated  electron 
disrrlbutions  depend  strongly  on  the  assumed  base  temperature 
but  do  not  depend  appreciably  on  the  small  temperature 
changes  assumed  along  a  line  of  force,  and,  that  to  a  first 
orde.=  ,  vertical  profiles  and  distributions  along  the  line 
of  force  computed  assuming  a  constant  temperature  are 
adequate . 
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The  ratios  (subsequently  referred  to  as  R)  of  the 
electron  densities  at  conjugate  points  In  the  two  hemispheres 
predicted  by  the  theory  were  calculated  for  a  series  of  heights 
(Fig.  15£tb).  The  temperature  distribution  along  a  line  of 
force  was  assumed  to  be  as  shown  In  Pig.  13t  the  temperature 
T  1  corresponding  to  that  at  the  base  level  In  the  northern 
(winter)  hemisphere.  If  it  is  supposed  that  the  compositions 
at  the  tase  level  are  the  same  in  winter  and  summer,  the 
ratio,  R  =  N(winter)/N(svmmer)  is  slightly  greater  than  unity 
and  tends  to  Increase  with  Increasing  latitude  (Fig. 15a).  If> 
Instead,  the  ionic  composition  at  500  Ion  is  assumed  to  depend 
on  the  temperature  as  given  by  Bauer  (1963)  (Pig.  8  and 
Equation  (65))  then  the  cuirves  of  Pig.  15b  are  obtained  and 
the  ratio  R  is  still  about  unity  except  for  the  1000  km 
level  for  which  it  decreases  from  1  at  the  equator  to 
approximately  0.5  at  high  latitudes. 

Fig.  16a, b,c  shows  N(h)  profiles  computed  assuming 
base  composition  1  for  the  two  assumed  temperature  models 
(see  Pig.  13  and  Equation  (66)).  It  was  assumed  that  the 
electron  density  at  1000  km  was  10^/cc.  Pig.  _l6a  corresponds 
to  an  equatorial  profile  (9  =  0°),  and  Figs.  lob.  c 
correspond  to  9  =  30°  and  60*^  respectively.  The  differences 
In  the  profiles  shown  In  Figs.  l6a,b,c  arise  only  because 
of  the  effect  of  the  centrifugal  force. 

In  order  to  convert  the  relative  electron  density 
distributions  to  absolute  values,  electron  density  data  for 
1000  km  obtained  from  the  Alouette  satellite  were  used 
(Thomas  and  Sader  1963).  Fig.  17,  based  on  the  results  of 
Thomas  and  Sader  (1963)  shows  the  mean  variation  of  the 
electron  density  at  1000  km  for  a  series  of  magnetically 
quiet  days  and  nights  In  summer  and  winter  as  a  function 
of  dip  latitude.  The  broken  lines  are  extrapolated  values. 

N(h)  profiles  were  computed  for  a  variety  of  conditions 
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using  the  data  of  Pig.  17  at  1000  km.  These  profiles  are 
shown  In  Figs.  l8a  -  d.  If  the  assumptions  about  temperature 
and  base  compositions  indicated  are  correct,  thetl  the  theorjf 
outlined  in  this  report  predicts  that  thd  electron  profiles 
in  summer  (day  and  night)  and  winter  (day  and  night)  are  aS 
shown  in  Figs.  l8a  -  d,  respectively. 

The  differences  between  the  curf/es  shown  in  Pig.  l8a  -  d 
arise  firstly  because  of  the  different  assumptions  about 
temperature  and  corapositlon^  and  secondly  because  of  changes 
in  the  electron  density  at  1000  knr  from  surnmar  to  winter 
and  from  day  to  night  used  in  converting  from  relative- 
distributions  to  absolute  distributions.  In  particular, 
the  rapid  fall  off  of  electron  density  with  height  in  the 
si-mmer  night  curves  is  mainly  attributable  to  the  rapid 
decrease  in  the  electron  density  at  1000  km  with,  dip 
latitude  over  the  range  50®  -  65° .  Typical  diurnal  and 
seasonal  changes  depicted  in  these  curves  are  illustra4;ed 
for  a  given  temperature  and  base  composition  In  Pig..  l3e 
which  shows  also  the  value  of  the  latitude  9^  corresponding' 
to  the  foot  of  the  field  line  (at  1000  km  above  the  eart^h) 
which  extends  to  a  given  distance  from  the  center  csf  tSiie 
earth  in  the  equatorial  plane. 

A  number  of  theoretical  N(h)  distributions  have  b€^n 
predicted  for  the  exospheric  plasma  and  these  are  dlscus:se'd 
in  Section  IV. C.  Some  typical  examples  are  reproduced  In 
Fig.  19a  (Bates  and  Patterson  I96I;  Dungey  195^;  Johnsoiri 
and  these  may  be  compared  with  the  curves  predicted  by  -the 
theory  presented  in  this  report. 

No  detailed  comparison  of  the  predicted  N(h)  profiles 
wltn  experimentally  observed  data  for  any  given  day  is 
possible  at  the  present  time.  However,  a  numoer  of  whistler 
mean  exospheric  electron  density  profiles  have  been  published 
covering  a  wide  range  of  observational  conditions,  A  number' 
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of  these  are  reproduced  In  Pig.  19^  and  are  labelled  a  to  h 
In  that  diagram  (Storey  1953#  Allcock  1959;  Carpenter  and 
Angerami  (quoted  by  Carpenter  1963);  Smith  'i960,  I96I; 

Pope  196! »  1962 j  Sohmelovsky  196O;  Schoute-Vanneck  and 
Muir  1963).  It  should  be  noted  that  the  observation  of 
Storey  (1953)#  qtidted  above,  refers  to  a  single  point 
only.  The  Incoherent  scatter  observations  of  Bowles  (1962) 
lie  within  the  arrow  indicated  in  Pig.  X9b. 

In  Pig:.  20  the  experimental  data  of  Pig.  19b  are 
superposed'  on  the  theoretical  curves  of  Pig.  19a.  The 
experiiffental  curves  of  Pig*  19b  cover  a  wide  variety  of 
observed  conditions  and  include  results  for  different  times 
of  day,  different  seasons  and  different  solar  epochs.  In 
general,  howeveiv  all  these  results  lie  within  the  approxi¬ 
mate  region  indicated  by  the  diagonally  shaded  area  in 
Pig.  ^0.  The  values  of  electron  density  measured  by 
Al;ouette  lie  within  the  approximate  zone  indicated  by  the 
vertically  shaded  area. 

It  is  clear  from  Fig.  20  that  when  the  values  of 
etectron  density  given  by  the  Alouette  data  for  summer 
nights  are  used,  the  theory  outlined  in  this  report  gives 
reasonably  good  agreement  with  the  observations  (curves  c  and 
whereas  the  curves  a,  f,  e,  b,  g  indicate  too  small  a 
(^crease  of  electron  density  with  height.  Furthermore  the 
present  theory  provides  a  betcer  agreement  with  the  observed 
sHope  of  the  N(h)  profile  near  the  level  of  the  Alouette 
orbit-.  It  should  be  noted  that  the  curve  labelled  e  represents 
probable  extreme  value  indicated  by  the  present  work  since 
it  corresponds  to  summer  daytime  conditions  when  the  assximp- 
tion  that  there  is  no  production  of  ionization  near  the  base 
level  is  most  likely  to  be  Invalid. 
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Until  further  information  is  available  concerning  the 
exospheric  temperature  and  composition,  it  is  probably 
unfruitful  to  consider  detailed  comparison  of  theory  with 
experiment.  However,  the  number  of  reasonable  choices  is  not 
unlimited  and  for  example  (Pig.  21),  the  N(h)  curve  of 
Smith  (i960,  1961)  is  well  matched  by  the  present  theory  if 
the  temperature  is  assumed  to  be  1000®K  and  the  starting 
electron  density  at  1000  km  Is  assumed  to  be  roughly  that 
corresponding  to  the  curve  labelled  "summer  night"  in 
Pig.  17.  Thus,  quite  good  agreement  can  bs  obtained  with 
values  which  might  be  reasonally  expected  to  apply  at  the 
time  when  the  whistler  observations  were  made. 

The  general  agreement  between  theory  and  experiment  Is 
illustrated  in  another  way  in  Pig.  22  which  shows  the- theoret-i^ 
ically  predicted  values  of  the  ratio  R  as  a  function  of 
geomagnetic  latitude,  together  with  observations  of  R  at 
1000  km  measured  by  the  Alouette  satellite.  The  shaded 
area  corresponds  to  values  within  which  all  the  observed 
values  of  the  ratio  R  lie  and  thus  allows  for  changes  from 
night  to  night,  etc.  The  circled  points  give  the  values 
for  magnetically  quiet  days. 


-  41 


C.EL-63-IIO 


r»rvMPaRTcinM  wtt’w  nmwpR  unRif  awn  stimmarv  nw  RRQTTT.mct 

v-v*  «-«*  %  WAV  **  <«»^  «*  w*  •*«•««>  «  v«  w**  V  w 


A.  DISCUSSION 

The  basic  concepts  which  govern  the  distributions  of 
ions  and  electrons  in  the  exosphere  have  been  discussed 
in  a  number  of  important  papers  in  the  scientific  literature. 
It  is  generally  agreed  that  the  main  constituent  of  the 
atmosphere  at  great  heights  above  the  earth  (above  approxi¬ 
mately  1500  kilometers)  is  hydrogen.  Hydrogen  ions  are 
formed  by  the  reaction 

0"^  +  H  +  0 


as  was  f.irst  pointed  out  by  Dungey.  The  relative  equilibrium 
concentrations  in  the  exosphere  are  determined  largely  by 
the  importance  of  this  reaction  near  the  base  of  the  neutral 
particle  exosphere  where  the  cross-section  for  collisions 
between  ions  and  neutral  particles  is  relatively  high.  The 
base  of  the  exosphere  at  approximately  550  kilometers  is 
believed  to  be  the  most  probable  level  of  origin  for  protons 
entering  the  exosphere  from  bolow  after  formation  in  the 
charge  exchange  reaction  referred  to  above  (Johnson  i960). 
Thus,  the  relative  abundances  of  O"^  and  (and  also  of  the 
Be"*"  ions)  at  the  base  of  the  exosphere  controls  the  electron 
and  ion  distributions  throughout  the  exosphere.  The  charge 
exchange  reaction  referred  to  above  proceeds  vei-^  rapidly 
near  the  peak  of  the  P2  layer  and  the  proton  number  density 
is  given  by  the  chemical  equilibrium  expression  (Hanson  and 
Ortehburger  196I) 


[R' 


m 


iiafiere  the  square  brackets  denote  the  concentrations  of  the 
constituents  bracketed. 
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It  has  also  been  shown  by  Hans  an  and  Crtenburger  (I961) 
that  the  oxygen  ions  near  the  F2  peak  are  only  weakly  coupled 
with  the  protons  in  the  protonosphere  and  it  seems  likely 
that  the  time  constant  for  the  establishment  of  equilibrium 
between  the  two  hemispheres  near  the  P  region  peak  is  greater 
than  one  day.  The  results  of  Hanson  and  Patterson  (1963) 
and  of  Hanson,  et  al  (1963)  provide  an  estimate  of  the  proton 
fluxes  likely  to  be  involved  as  the  result  of  the  diurnal 
variation  of  the  abundance  of  the  atomic  hydrogen  in  the 
exosphere.  They  considered  the  flow  of  hydrogen  into  and  out 
of  the  protonosphere  by  the  charge  exchange  process  referred 
to  above,  the  escape  of  hydrogen  from  the  daytime  exosphere 
associated  with  hi^er  daytime  temperatures  and  also  the 
lateral  flow  of  hydrogen  around  the  earth  due  to  any  asymmetry 
in  the  distribution  of  hydrogen  around  the  earth. 

It  is  clear  from  this  work  and  from  other  considerations 
that  some  modifications  to  a  simple  equilibrii”^  theory  are 
probably  required  to  describe  completely  the  physics  of  the 
exosphere  in  terms  of  diurnal  and  other  changes.  However, 
before  a  non-equilibrium  theory  is  properly  applied,  it  is 
valuable  to  explore  the  extent  to  which  an  equilibrium  theory 
is  adequate  to  explain  the  observational  data  and  this  has 
been  the  purpose  of  the  work  reported  herein.  In  particular, 
it  seems  likely  that  an  equilibrium  theory  might  be  applicable 
for  explaining  seasonal  variations  which  are  likely  to  occur 
in  a  given  hemisphere  as  well  as  other  changes  involving  time 
constants  greater  than  one  day. 

The  structure  of  the  lower  exosphere  has  been  discussed 
in  some  detail  by  Hanson  (1962),  Bauer  (I962,  I963)  and  by 
Gliddon  (1962)  who  considered  the  distributions  of  the  0^, 

He"^,  and  H^  ions  together  with  the  electrons  under  diffusive 
equilibrium.  The  corresponding  distributions  at  much  greater 
altitudes  have  been  discussed  by  Dungey  (1954),  Johnson  (i960) 
and  by  Bates  and  Patterson  (I961), 
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The  main  differences  between  the  work  presented  in  this 
report  and  the  earlier  work  referred  to  above  relate  to  the 
following  points; 

1.  Under  the  circumstances  where  diffusive  equilibrium 
applies,  it  is  well  known  that  whereas  in  a  non-ionized 
atmosphere  each  constituent  is  distributed  independently 

of  the  other  constituents,  the  distribution  of  a  light 
ionic  constituent  in  diffusive  equilibrium  ^  dependent 
upon  the  presence  of  the  others  because  of  the  electric 
field,  E,  which  arises  from  the  slight  charge  separation 
between  electrons  and  positive  ions  (Dungey  195^^,  flange  i960, 
Hanson  and  Ortenburger  1961).  In  the  present  work,  it  has 
been  assumed  that  the  electrons  and  ions  can  diffuse  only 
along  the  geomagnetic  field  lines  and,  therefore,  the 
direction  of  the  electric  field,  E,  has  been  assumed  to  be 
along  the  particular  field  line  concerned. 

2.  It  is  shown  that  it  is  convenient  to  work  in  terms 

of  a  parameter,  z,  called  the  temperature -modified  geopotential 
height.  The  equations  have  been  developed  to  allow  for 
the  possibility  that  the  electron  and  ion  temperatures  may 
be  different  and  that  a  temperature  gradient  may  exist  along 
a  line  of  force. 

3.  As  was  first  pointed  out  by  Bates  and  Patterson 

(1961),  up  to  the  present  time,  the  theories  of  electron 
and  ion  distributions  in  the  earth's  exosphere  have  not 
taken  into  account  conditions  in  the  upper  part  of  the 
earth's  ionosphere  which  may  be  strongly  dependent  on 
latitude.  Such  a  latitude  dependence  would  influence  the 
electron  density  distribution  in  the  equatorial  plane  since 
the  diffusion  of  charged  particles  in  the  earth's  exosphere 
is  appreciable  only  along  geomagnetic  field  lines.  Since 
our  N(h)  profiles  have  been  deduced  from  the  theoretically 
computed  N(s)  distributions  (i.  e.,  the  calculated  equili¬ 
brium  distributions  along  a  line  of  force  when  diffusion 
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along  the  field  lines  is  operative),  the  starting  electron 
density  assvmed  at  the  base  level  and,  in  particular,  its 
variation  with  latitude,  is  an  extremely  important  factor. 

An  examination  of  the  results  of  Thomas  and  Sader  (1963) 
shows  that  the  electron  density  at  1000  km  as  observed  by  the 
Aloiette  satellite  does  vary  quite  strongly  with  latitude  and, 
in  particular,  there  is  an  important  minimum  in  the  curve  under 
certain  circumstances  near  a  dip  latitude  of  65°.  This  latitude 
dependence  was  inserted  into  the  theory  so  that  the  actual 
observed  electron  densities  at  1000  km  were  used.  The  computed 
N(h)  curves  in  the  equatorial  plane  at  great  distances  from 
the  earth  were  then  compared  with  observations  from  whistler 
data.  It  was  found  that  the  disagreement  between  theory  and 
observations  of  the  exospheric  plasma  electron  density  distri¬ 
bution  in  the  vicinity  of  the  Alouette  orbit  (see  Pig.  20) 
was  removed.  Under  certain  circumstances,  the  disagreement  at 
greater  distances  from  the  earth  was  also  removed  when 
reasonable  assumptions  about  tne  exospheric  temperature  and 
about  the  relative  ion  densities  at  500  krn  were  made  and  the 
Alouette  data  for  summer  nights  used.  It  is  possible  that 
at  the  times  when  the  N(h)  slopes  far  out  do  not  agree  with 
theory,  that  some  of  the  assumptions  made  in  the  theory  were 
invalid.  In  particular,  it  might  be  the  case  that  composition 
and  temperature  changes  at  1000  km  were  inapplicable,  or  that 
the  level  above  which  diffusive  equilibrium  occurs  varies  with 
time. 

The  slope  of  the  N(h)  profile  in  the  lower  exosphere  (approxi¬ 
mately  600  -  1000  km). 

It  has  been  assumed  in  the  theory  that  the  production 
of  electrons  by  the  sun's  ionizing  radiations  and  the  loss 
of  electrons  by  recombination  are  both  non-existent.  Thus 
no  attempt  is  made  to  discuss  the  N(h)  distributions  at  or 
near  the  peak  of  the  P2  layer,  though  a  great  deal  of 
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observational  and  theoretical  data  are  available  (see  for 
example  Thomas  1963,  Rishbeth,  Lyon  and  Peart  1963  and 

i.  ^  i.  1.  V  t  I.  ^  xt  t  /  •  \m&x  o  v^o  o  X  VI :  Aico 

been  confined  to  consideration  of  the  theory  which  might 
apply  above  the  critical  level  (Johnson  i960)  at  about  500 
km  abc ve  the  earth’s  surface. 

Bauer  (1962),  (1963)>  assumed  that  in  the  vicinity  of 

the  Alouette  orbit  there  is  diffusive  equilibrium  with  O"^, 

+  + 

He  and  H  as  the  main  constituents  and  later  was  able  to 
fit  a  topside  profile  from  the  Alouette  satellite  on  the 
basis  of  an  equilibrium  theory  down  to  quite  low  altitudes. 

The  results  of  the  theory  outlined  in  the  present  work  in 
the  vicinity  of  the  Alouette  orbit  agree  with  Bauer’s  slopes. 

The  main  differences  between  the  work  presented  in  the 
report  and  that  referred  to  above  is  that  Bauer  considered 
an  Isothermal  case  and  does  not  allow  for  centrifugal  force: 
both  assumptions  being  Justified  for  the  cases  which  Bauer 
considers  since  he  restricts  his  theory  to  distances  which 
do  not  exceed  about  one  earth  radii  above  the  P2  region 
peak.  Bauer,  however,  considers  equilibrium  in  a  vertical 
column  rather  than  in  a  column  confined  by  a  tube  of  force. 

In  the  approach  outlined  in  this  report,  a  vertical  profile 
is  very  strongly  influenced  by  a  north-south  gradient  in  the 
temperature  and  ionic  composition  at  the  base  level  or  by  a 
north-south  gradient  in  the  electron  density  at  the  base  level. 

Since  the  calculations  described  in  this  report  were 
carried  out,  a  number  of  N(h)  profiles  deduced  from  the 
Alouette  satellite  data  have  been  circulated.*  The  available 

*  ftadlo  Research  Station,  Slough,  England,  Topside  N(h) 

Profiles,  Vol.  I,  Dec.  1963;  King,^J.  W.,  et  al.  Radio 

Research  Report  No,  IM  9k,  July  1963;  and  King,  J.  W,, 

et  al.  Report  No.  IM  112,  Dec.  I963. 

* 
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profiles  correspond  to  equatorial  crossings  of  the  satellite 
and  were  observed  at  the  telemetry  station  Singapore.  Prom 
the  N(z)  plots  of  King,  et  al  (1963)  (see  footnote,  pg.  46), 
measurements  were  made  of  the  scale  height  for  electrons  at  1000 
km  for  one  daytime  and  one  nighttime  revolution:  26  November 
1962,  14:58  local  time,  and  I5  April  1963,  20:28  local  time, 
respectively. 

For  the  daytime  pass,  the  scale  height  at  1000  km  varies 
between  170  km  at  14®N  latitude  (consistent  with  O"^  .at  1380°K) 
and  260  km  at  12°S  latitude  (consistent  with  a  mixture  of  O"^ 
and  He"^  at  1380°K).  Around  the  latitude  l4°N  the  electron 
distribution  is  a  pure  exponential  down  to  350  km,  supporting 
the  idea  that  0"^  is  indeed  strongly  predominant. 

For  the  nighttime  pass,  the  scale  height  at  1000  km  varies 
between  38O  km  at  20°N  (consistent  with  strongly  predominant 
He"^  at  780°K--or  a  mixture  of  the  He'^’with  a  significant 
amount  of  0*^,at  a  slightly  higher  temperature)  and  620  km  at 
9°S  (consistent  with  strongly  predominant  He"^  at  1270°K--or 
a  mixture  of  He^  and  a  significant  amount  of  at  a  slightly 
lower  temperature).  The  preceding  measurements  show  then 
that  the  lower  boundary  of  the  helium  layer  is  lower  at 
night  than  during  the  day. 

Prom  the  above  it  appears  that  the  scale  height  for 
electrons  at  1000  km  varies  strongly  from  day  to  night  (by 
a  factor  2.4  in  the  above  measurements)  and  also  varies 
with  latitude.  It  is  apparent  that,  at  1000  km,  both  the 
ionic  composition  and  the  temperature  must  vary  diurnally 
and  also  with  latitude,  at  least  in  the  equatorial  regions. 

Also,  it  is  important  to  note  that  the  assumption  that  H'*’ 
strongly  predominates  at  1000  km  would  lead  to  extremely 
low  temperatures  in  order  to  explain  scale  heights  as  low 
as  170  km  (T  less  than  100°K). 
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The  slope  of  the  profile  In  the  upper  exosphere . 
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shows  that  at  vevy  gi^eat  heights  the  rate  at  which  the 
electron  density  decreases  with  height  is  often  less  than 
the  experimental  observations  indicate.  Furthermore,  the 
slopes  of  the  profiles  near  1000  km  for  curves  a,  b,  f,  g. 
Pig.  19a f  are  not  in  agreement  with  those  observed  from 
electron  density  profiles  deduced  in  the  vicinity  of  the 
Alouette  orbit.  The  theoretical  profiles  of  Dungey  {195^) > 
Bates  and  Patterson  (196I)  and  Johnson  (i960)  allow  onlj'^ 
for  the  presence  of  and  not  for  helium.  Thus,  the  slopes 
are  small  even  at  low  altitudes.  The  present  work  produces 
slopes  in  agreement  with  those  observed  by  the  Alouette 
satellite  near  its  orbit  in  conformity  with  the  results  of 
Bauer. 


At  the  present  time,  until  more  information  is  available 
about  the  positive  ion  abundances  at  1000  km,  It  is  not  known 
to  what  extent  profiles  observed  at  all  times  can  be  fitted 
by  a  simple  equilibrium  theory.  However,  it  is  clear  that 
the  main  features  of  the  seasonal  variations  will  probably, 
to  a  large  extent,  be  governed  by  the  broad  precepts  dictated 
by  an  equilibrium  theory,  although  there  will  be  diurnal 
changes  superimposed  on  the  broad  seasonal  changes.  The 
theory  outlined  in  the  present  paper  does,  under  certain 
circumstances,  lead  to  electron  density  profiles  which  fail 
off  at  great  altitudes  at  a  rate  which  is  in  agreement  with 
the  experimental  data  obtained  from  whistlers  and  at  the 
same  time  gives  the  correct  slope  near  the  Alouette  orbit. 
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B,  CONCLUSIONS 


1.  The  equilibrium  distributions  of  electrons  and  ions 
in  the  earth's  exosphere  have  been  computed  for  a  wide  variety 
of  conditions  on  the  assumption  that  diffusive  equilibrium 
applies . 

2.  The  theory  assumes  that  diffusion  of  electrons  and 
ions  is  possible  only  along  the  geomagnetic  field  lines. 

3.  The  presence  of  He”^  ions  near  the  base  of  the 
exosphere  has  been  allowed  for  in  the  calculations.  This  is 
vital  for  an  understanding  of  the  changes  accompanying  the 
transition  from  day  to  night  conditions. 

4.  The  most  crucial  factors  governing  the  form  of  these 
distributions  are  the  exospheric  temperature  and  the  ionic 
composition  at  the  base  level  assumed  to  be  at  500  kilometers 
above  the  earth's  surface. 

5.  The  strong  latitude  dependence  of  electron  density 
at  1000  km  as  observed  by  the  Alouette  satellite  has  been 
taken  into  consideration  in  the  calculations.  As  a  result,  it 
Is  clear  that  the  discrepancy  previously  observed  between 
theory  and  experiment,  namely,  that  the  electron  density 
gradient  at  great  heights  predicted  by  the  theories  turned 
out  to  be  much  smaller  than  that  indicated  by  the  whistler 
observations,  is  now  removed  under  certain  circumstances. 

6.  The  effect  of  the  earth’s  centrifugal  force  is 
included  in  the  calculations — its  neglect  would,  at  the 
high  altitudes  lead  to  electron  density  values  which  are  too 
low  at  a  given  height  by  an  amount  which  does  not  exceed 

20  percent,  for  altitudes  less  than  about  7  earth  radii. 

Although  it  is  not  suggested  that  an  equilibrium  theory 
is  adequate  to  describe  the  entire  properties  of  the  plasma 
distribution  in  the  exosphere,  it  is  quite  clear  that  a  first 
order  agreement  can  be  obtained,  at  any  rate,  in  terms  of 
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seasonal  dependencies .  The  basic  questlc 
density  at  the  Alouette  orbit  varies  with 
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Further 


n  of  why  the  electron 
latitude  in  the  way 
work  Is  needed  to 


determine  the  nature  and  the  relative  abundances  of  the 
positive  ions  at  the  base  of  the  exosphere  and  the  dependence 
of  these  relative  abundances  on  the  exospheric  temperature. 
Measurements  of  the  relative  abundances  on  a  dally  routine 
basis  through  satellite  observations  will  probably  be  necessary 
before  it  is  possible  to  determine  accurately  the  circumstances 
in  which  an  equilibrium  theory  becomes  inapplicable. 
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APPENDIX  A,  GEOMETRY  OP  THE  DIPOLE  FIELD 


It  was  assumed  in  this  report  that  the  actual  magnetic 
field  of  the  earth  could  be  approximated  by  a  centered 
dipole  field,  as  indicated  in  Mlodnosky  and  Helliwell  (1962). 
Most  of  the  equations  quoted  in  this  appendix  appear  in 
Pig.  1  of  that  paper,  to  which  reference  should  be  made. 

Thus,  the  equation  of  a  line  of  force  in  a  dipole  field  is 


r  cos^e 
o  COS  8 


(Al) 


in  which  the  quantities  are  as  shown  in  Pig.  1,  page  4,  of 
this  report.  In  this  figure,  the  distance  s^  along  the  line 
of  force  with  feet  at  G^,  from  the  equatorial  plane  to  the 
latitude  e  is 


Si(e^,0) 


_ ^o  (x  +  sinhx.coshx) 

2  yPS  cos^e^ 

°  (A2) 

- — 5 —  lln  (y  +  VI  +  y^)  +  y  y/l  +  y^] 

2  ys  cos^e^ 


where 


y  =  sinhx  =  yj  sinG 


We  called  Sjyj(6^)  the  value  of  s^^  corresponding  to 


=  G  ,  as  shown  in  Pig.  7,  page  29: 

w 


S„(9„)  = 


2  JJ 


COS  G 


fin  (y^+  >/i  +  y^)+  y^i  +  y^J  (A3) 


9 


with 


y,3  -  yi  Eine^ 


For  con.enience,  the  distance  s(9^jG)  botv/een  a  bare 
level  (e^)  and  a  latitude  0  was  used  (see  Pig.  1,  page  i. ) : 

=  (a2*) 


The  increment  of  distance  ds,  corresponding  to  an 
increment  de  (keeping  6^  constant)  is  easily  shown  to  be 


ds 


r 

o 


2 

cos  0 

o 


3  sin^e 


cose  d9 


It  is  clear  that  by  means  of  eq.  (Al)  the  distances 
^  be  written  as  functions  of  (0^,r),  and 
this  is  convenient  for  electron  and  ion  density  calculation;, 
in  which  the  ^-xospheric  temperature  is  assumed  constant. 

This  procedure  was  used  in  the  ccmputatlons  leading  to 
Figs.  12  and  l8 . 

It  should  be  noted  chat  the  inversion  of  eq.  (42) 

[or  (A4)]  to  give  9  as  a  function  of  9^  and  s^  (or  s)  is 
only  possible  numerically.  Thus,  in  the  computations  leading 
to  Pi£;s.  1^  -  i6,  in  which  the  temperature  is  assumed  a 
function  of  the  distance  s,  the  numerical  integrations  t 
get  z  (from  eq.(50))  were  made  using  9  as  the  independent 
variable . 
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APPENDIX  B.  NOTE  ON  THE  ION  DISTRIBUTION  CURVES 

For  an  isothermal  exosphere,  the  relative  electron  and 
ion  densities  can  be  conveniently  plotted  as  functions  of  z. 
Assuming  that  the  electrons  and  ions  all  have  the  same  tempera¬ 
ture  (i.  e.,  C  =  1),  their  relative  densities  ai?e- obtained 
from  equations  (31)  -  (33)  respectively  (pages  l6  -  17): 


ng(z) 


exp  (-z/Hj^)  +  Hg  exp  (-z/Hg)  +  exp  (-z/Hj)!? 


(Bl) 


nj^(z) 


(B2) 


’’l _ ^  (  z  1  1  '’eo  ) 


(B3) 


If  one  ion  is  strongly  predominant,  so  that  only  the 
term  for  this  particular  ion  is  important  on  the  right-hand 
side  of  equation  (Bl),  then. 


■I”  -  Vr  “p  I'  Iq;) 


(Bit) 


where  the  suffice  i  refers  to  the  strongly  predominant  ion. 
Equation  (b4)  is  a  simple  exponential  and  is  represented  by 
a  straight  line  on  serai-logarithmic  paper,  with  slope  2H^, 
(twice  the  scale  height  of  the  non-ionized  specimen).  These 
straight  lines  are  drawn  in  Pigs.  9a  -  d  for  O'*',  He**", 

(c.f.  Hanson  (1962))  and  are  labelled  2H^,  2H2  and  BH^ 
respectively. 
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Reference  to  the  N(z)  curve  of  Fig.  gb  shows  that  the 
simple  exponential  behaviour  of  electrbn  density  referred  to 
above  is  exhibited  below  100  km  where  O"^  is  strongly 
predominant  (note  that  the  ordinate  is  z  riot  h)  and  is  also 

4. 

exhibited  above  I5OO  km  where  H  is  strongly  predominant. 
Around  500  km,  although  He"*"  is  predominant,  it  is  not 
sufficiently  strongly  predominant  and  a  straight  line 
olstributlcn  is  not  achieved. 

Comparison  of  Figs.  9c  and  d  shows  the  effect  in  the 
slope  of  the  N(z)  curve  of  changing  the  composition  at  the 
base  level.  According  to  Table  1,  Fig.  9d  corresponds  to  a 
smaller  percentage  of  He"^  and  H"'’  at  the  base  level,  so  that 
the  o'*"  is  predominant  over  a  greater  height  range.  Thus, 
the  electron  density  maintains  a  high  rate  of  decrease  over 
a  greater  height  range  so  that  the  electron  density  turns  out 
to  be  smaller  for  the  same  z,  in  curve  d  than  in  curve  c 
(Fig.  10).  As  a  result,  one  sees  that  a  change  in  composi¬ 
tion  is  for  more  important  than  a  change  in  temperature  in 
determining  the  value  of  the  relative  electron  density  as 
a  function  of  z  (Fig.  10). 

In  equation  (B2),  the  term  exp  ( -z/H^ )  decreases  with 
increasing  z,  whereas  the  term  n^^/n^{z)  increases  with 
Increasing  z.  It  is  therefore  possible  for  the  ion  density 
to  Increase  with  z,  since  it  is  made  up  of  the  product  of 
these  two  factors. 

The  ratio  of  densities  of  two  different  ions  as  function 
of  z  may  be  obtained  from  equation  (B3): 


n^(z) 

FTjTTr  “ 


(B5) 


which,  in  the  logarithmic  plot  of  Figs.  9a  -  e,  is  a  straight 
line  with  (positive)  slope  Hj  -  Hj)  passing  through  the 
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0 


point  at  z  =  0.  These  lines  are  not  shown  in  Pigs,  ya  - 

but  may  readily  be  Included. 

Equation  (B3)  may  be  written  in  a  more  convenient  form: 


n^(z) 

^i 

n 

n 

eo 

eo 

(B6) 


so  that  the  geometric  mean  of  the  electron  and  1th  ion 
densities  relative  to  the  electron  density  at  the  base  level 
is 


which  is  o  ir  previous  equation  (B4,)  giving  the  asymptotic 
behaviour  for  the  relative  electron  density  curve  when  the 
ith  ion  is  strongly  predominant.  Thus  from  (b6) 


The  left-hand  side  of  this  equation  is  the-  arithmetic  mean 
of  the  logarithms  of  the  electron  density  and  the  density  of 
the  ith  ion  relative  to  the  electron  density  at  the  base 
level. 

Thus,  if  relative  electron  density  curves  are  plotted  using 
a  logarithmic  scale  in  density  (as  in  Figs.  9^  -  e )  then  the 
distribution  curve  for  the  ion  i  may  be  derived  from  a  simple 
geometric  construction  as  follows.  First  the  straight  line 
(B^)  is  plotted  (slope  ,  crossing  the  abscissa  at  a  value 
/"H )  then  the  relative  electron  distribution  (B1)  Is 
plotted.  For  each  value  of  z,  a  point  is  plotted  such  that 
the  horizontal  distance  from  the  point  to  the  straight  line 
is  the  same  as  that  from  the  straight  line  to  the  electron 
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distribution.  The  locus  of  these  points  is  the  ion 
distribution  curve  (Pig.  9b).  This  procedure  may  be 
repeated  for  each  ion. 

It  is  clear  from  this  construction  that  the  ionic 
transition  levels  (the  values  of  z  at  which  two  ions  have 
the  same  density,  i.  e.,  the  intersection  point  of  two 
ionic  distribution  curves)  may  be  alternatively  identified 
as  the  Intersection  of  the  corresponding  asymptotic  straight 
lines. 
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ASSUMPTIONS  ABOUT  COMPOSITION  AT  500  KM 

COMPOSITION  1 

COMPOSmOKS  2 

A 

[T  =  constant] 

1)2=  2.0  X  10‘^ 
1)^  =  1.6  X  10"^ 

T  “(K) 

’'2 

^5 

2a 

2b 

2c 

2d 

2e 

500 

T50 

1250 

1500 

2000 

6.5  X  10“^ 

5.6  X  10"^ 

1.1  X  10 

6.5  X  lO"^ 
2.0  X  10"^ 

1.6  X  lO’^ 

5.1  X  lO"^ 

-U 

5.1  X  10 

~k 

1.6  X  10 

1.6  X  lO'^ 

B 

[T  =  T(s)] 

-2 

Tig  =  2.0  X  10 

T^  =  1.6  X  10"^ 

2f 

Tl2 

T]5 

lliOO 

8.0  X  10*^ 

2.6  X  10 

TABLE  1 

TEMPERATURE  AND  COMPOSITION  ASSUMPTIONS 
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ASSUMPECONS  ABOUT  TEMPERATURE 


TABLE  2.  -SUMMARY  OF  PRESENTATION  OP  RESULTS 
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ASSUMPTIONS  ABOUT  COMPOSITION  AT  500  KILOMEl’ERS 
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RELATIVE  POSITIVE  ION  DENSITIES.  BAUER  1963 


z,  km 


RELATIVE  DENSITIES 


FIG.  9a.  RELATIVE  ELECTRON  AND  ION  DENSITIES  AS  A  FUNCTION  OF  THE 
GEOPOTENTIAL  HEIGHT,  z,  FOR  AN  ISOTHERMAL  EJ^OSPHERE,  The  form  of 
the  distributions  dej>9nds  strongly  on  the  assumed  base  level  com¬ 
position  and  the  assu^d  temperature  as  may  be  seen  by  cojnparing 
the  above  with  the  subsequent  figures. 
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RELATIVE  DENSITIES 


FIG.  9b.  RELATIVE  ELECTRON  AND  ION  De(SITIES  AS  A  FUNCTION  OF  THE 
GEOPOTENTIAL  HEIGHT,  z,  FOR  AN  ISOTHERMAL- EXOSPfflSRE.  Conparison 
with  Fig.  9c  indicates  the  effect  on  the  distributions  of  chang¬ 
ing  the  temperature  only. 
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BASE  COMPOSITION  I 


FIG.  9c.  RELATIVE  ELECTRON  AND  ION  DENSITIES  AS  A  FUNCTION  OF  THE 
GEOPOTENTIAL  HEIGHT,  z,  FOR  AN  ISOTHERMAL  EXOSPHERE.  Comparison 
with  Fig.  9b  indicates  the  effect  on  the  distributions  of  chang¬ 
ing  the  temperature  only.  Comparison  with  Fig.  9d  reveals  the 
effect  of  changing  the  base  composition  only. 
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RELATIVE  ELECTRON  DENSITY 


FIG.  10.  SUNMARY  OF  THE  RESULTS  OF  FIGS.  9a  >  e  FOR  THE  RELATIVE 
ELECTRON  DENSITY  DISTRIBUTIONS  IN  THE  EXOSPHERE  OVER  A  WIDE 
RANGE  OF  CONDITIONS. 
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ION  TRANSITION  LEVELS,  Z,  km 
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FIG.  11.  THE  VARIATION  WITH  TEMPERATURE  OF  THE  LEVELS  AT  WHICH 
THE  IONS  INDICATED  HAVE  EQUAL  ABUNDANCES.  The  broken  lines 
correspond  to  the  cases  when  the  composition  at  the  base  level 
depends  on  the  temperature  there  (through  equations  (65)).  The 
continuous  lines  correspond  to  the  case  when  the  composition  at 
the  base  level  is  independent  of  temperature  (Composition  l). 
These  results  are  in  agreement  with  those  presented  earlier  by 
Bauer  (l962). 
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500  KM  ABOVE  THE  EARTH)  AT  GEOMAGNETIC  LATITUDE  45®  (SEE  TABLES  1 
AND  2 )  . 
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FIG.  12b.  TO  ILLUSTRATE,  FOR  AN  ISOTHERMAL  EXOSPHERE,  THE  RELATIVE  ELECTRON 
DENSITY  DISTRIBUTIONS  AIX)NG  A  FIELD  LINE  WITH  FEET  (AT  500  KM  ABOVE  THE 
^TH)  at  geomagnetic  latitude  55“  (see  TABLES  1  AND  2).  The  broken  lines 
.ndicate  the  distributions  which  would  have  been  obtained  if  the  centrifugal 
force  had  been  neglected. 
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FIG.  13.  SCHEMATIC  DIAGRAM  TO  ILLUSTRATE  THE  VARIATION  OF  TEMPERATURE  ALONG  A  LI.NE  OF  FORCE. 
Ti  (s)  is  used  for  the  temperature  distributions  in  the  northern  (winter)  hemisphere  and 
T2(s)  for  the  southern  (summer)  hemisphere. 
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FIG.  14a.  TO  ILLUSTRATE.  FOR  A  NON- ISOTHERMAL  EXOSPHERE  (SEE 
FIG.  13).  THE  RELATIVE  ELECTRON  DENSITY  DISTRIBUTION  ALONG  A 
FIELD  LINE  WITH  FEET  (AT  500  KM  ABOVE  THE  EARTH)  AT  GEOMAG¬ 
NETIC  LATITUDE  35**  (SEE  TABLES  1  AND  2). 
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FIG.  14b.  TO  ILLUSTRATE,  FOR  A  NON- ISOTHERMAL  EXOSPHERE  (SEE  FIG.  13),  THE  RELATIVE 
ELECTRON  DENSITY  DISTRIBUTION  ALONG  A  FIELD  LINE  WITH  FEET  (AT  500  KM  ABOVE  THE 
EARTH)  AT  GEOMAGNETIC  LATITUDE  45°  (SEE  TABLES  1  AND  2). 
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FIG.  14c,  TO  ILLUSTRATE,  FOR  A  NON- ISOTHERMAL  EXOSPHERE  (SEE  FIG.  13),  THE 
RELATIVE  ELECTRON  DENSITY  DISTRIBUTION  ALONG  A  FIELD  LINE  WITH  FEET  (AT 
500  KM  ABOVE  THE  EARTH)  AT  GEOMAGNETIC  LATITUDE  55°  (SEE  TABLES  1  AND  2), 
The  broken  lines  indicate  the  distributions  which  would  have  been  obtained 
if  the  centrifugal  force  had  been  neglected. 
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FIG.  15b.  THE  PREDICTED  RATIO  OF  THE  ELECTRON  DENSITY  IN  WINTER  TO  THAT  AT  THE  SAME  LOCATION 
AND  ALTITUDE  IN  SUMMER  FOR  A  NUMBER  OF  DIFFERENT  ALTITUDES.  Compositions  1  and  2f  (See 

Table  1)  are  used  at  500  km,  as  was  done  previously.  The  assumed  temperature  distribution 
is  illustrated  in  Fig.  13. 
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FIG.  16a.  EQUATORIAL  N(h)  PROFILE.  The  electron  density  at  1000  km  was 
assumed  to  be  lO^/cc  (See  Tables  1  and  2).  The  temperature  distributions 
along  a  field  line  are  given  by  equation  (66)  and  are  illustrated  schemati¬ 
cally  in  Fig.  13. 
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FIG.  16b.  N(h)  PROFILE  AT  GEOMAGNETIC  LATITUDE  30®.  The  electron  density  at 
1000  km  was  assumed  to  be  104/cc  (See  Tables  1  and  2).  The  temperature 
distributions  along  a  field  line  are  given  by  equation  (66)  and  are  illus¬ 
trated  schematically  in  Fig.  13. 
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HEIGHT,  10^  km 

FIG.  16c.  N(h)  PROFILE  AT  GEOMAGNETIC  LATITUDE  60°  .  The  electron  density  at 
1000  km  was  assumed  to  be  lO^/cc  (See  Tables  1  and  2).  The  temperature  dis¬ 
tributions  along  a  field  line  are  given  by  equation  (66)  and  are  illustrated 
schematically  in  Fig.  13, 
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FIG,  I8a.  THEORETICALLY  PREDICTED  EQUATORIAL  N(h)  PROFILES  BASED  OX  OBSERVED 
AU)UETTE  DATA  FOR  THE  ELECTRON  DENSITY  AT  1000  KM  FOR  SUMMER  DAYS.  Curves  for 

wide  variety  of  isothermal  exospheric  temperatures  and  compositions  are  illus¬ 
trated  (See  Tables  1  and  2l. 
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FIG.  18b.  THEORETICALLY  PREDICTED  EQUATORIAL  N(h)  PROFILES  BASED  ON  OBSERVED 
ALOUETTE  DATA  FOR  THE  ELECTRON  DENSITY  AT  1000  KM  FOR  SUMMER  NIGHTS.  Curves  for 
a  wide  variety  of  isothermal  exospheric  temperatures  and  compositions  are  illus¬ 
trated  (See  Tables  1  and  2).  Note  the  marked  increase  in  the  rate  of  fall  off 
of  electron  density  with  height  at  great  distances  from  the  earth. 
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FIG.  18d.  THEORETICALLY  PREDICTED  EQUATORIAL  N(h)  PROFILES  BASED  ON  OBSERVED 
ALOUETTE  DATA  FOR  THE  ELECTRON  DENSITY  AT  1000  KAI  FOR  WINTER  NIGHTS.  Curves  for 
a  wide  variety  of  isothermal  exospheric  temperatures  and  compositions  are  illus¬ 
trated  (See  Tables  1  and  2). 
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EXPERIMENTAL  PROFILES 

a  STOREY 


G.  19b.  EX/VMPLES  OF  EXPERIMENTAL  OBSERVATIONS  OF  EXOSPHERIC  ELECTRON  DENSITY 
( equator  I a L  prof I LES ) . 
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FIG.  20.  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  EXOSPHERIC  N(h)  PROFILES  IN  THE 
EQUATORIAL  PLANE.  The  curves  a  -  g  are  those  described  In  Fig.  19a.  The  diagonally 
shaded  area  defines  approximately  the  region  within  which  the  experimental  profiles 
of  Fig.  19b  are  observed  to  lie.  The  vertical  .shading  near  1000  km  indicates  the 
approximate  region  in  which  the  Alouette  observations  for  1000  km  lie. 
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G.  21.  COMPARISON  OF  THEORY  AND  EXPERIMENT.  The  theoretical  curve  was  computed 
assuming  a  constant  exospheric  temperature  of  1000  "k  and  Ionic  Composition  1 
VTable  l)  at  the  base  level.  The  electron  density  at  1000  km  was  taken  to  be 
that  given  by  Alouctte  for  summer  night  conditions  (Fig.  17), 
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FIG.  22.  OBSERVED  AND  PREDICTED  RATIOS  OF  THE  ELECTRON  DENSITY  .\T  1000  KM  AT  CON¬ 
JUGATE  POINTS  IN  THE  WINTER  AND  SUMMER  HEMISPHERES.  The  middle  line  giving  the 
observed  average  quiet  day  ratios  is  based  on  summer  and  winter  Alouette  observa 
tlons  at  Stanford.  The  ratio  is  approximately  0.5  over  the  latitude  range  con¬ 
sidered.  The  theoretically  predicted  ratios  are  in  reasonably  good  agreement. 


